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ABSTRACT 


The  dynamical  laws  governing  spinning  multipole  test 
particles  and  polarized  media  with  internal  spin  are  derived 
from  both  variational  principles  and  the  multipole  formalism 
of  extended  bodies. 

The  general  form  of  the  Lagrangian  equations  of  motion 
is  derived  for  a  spinning  multipole  particle  in  given 
external  fields.  We  then  consider  the  dynamics  of  a  contin¬ 
uous  medium  with  internal  spin  and  multipole  structure.  From 
a  four -dimensional  action  integral  the  field  equations 
relating  the  fields  generated  by  the  medium  to  its  bulk 
properties  are  derived,  together  with  the  balance  laws 
expressing  conservation  of  total  four -momentum  and  spin. 

A  natural  splitting  of  the  total  energy -momentum  tensor 
into  matter  and  field  parts  is  adopted  that  leads  to  a 
"generalized  Minkowski"  electromagnetic  energy  tensor.  In 
both  the  electromagnetic  and  the  gravitational  field  equa¬ 
tions  the  source  terms  contain  "polarization"  contributions. 

It  is  shown  that  the  multi  pole  formalism  may  be  used  to 
formulate  the  same  equations  of  motion,  balance  laws  and 
decomposition  of  total  energy -momentum  as  those  resulting 
from  variational  principles. 
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1.  INTRODUCTION 


This  thesis  formulates  the  dynamics  of  interacting  material 
and  fields,  where  the  material  possesses  inner  structure  in 
the  form  of  spin  angular  momentum,  electric  and  magnetic 
dipole  moments,  gravitational  quadrupole  moments  and  higher 
moments.  The  aim  is  to  derive  covariant,  general  relativ¬ 
istic  field  equations  coupling  the  gravitational,  electro¬ 
magnetic  and  other  fields  to  their  sources,  and  to  obtain 
balance  laws  expressing  conservation  of  total  charge, 
energy -momentum  and  angular  momentum.  The  major  task  will  be 
to  find  the  general  form  of  the  gravitational  source,  the 
total  energy -momentum  tensor  of  the  material  and  the  non- 
-gravi tat iona 1  fields. 

Prior  to  the  last  decade,  the  structure  of  the  energy- 
-momentum  tensor  was  mainly  of  academic  interest  to  general 
relativists.  Research  has  focused  mainly  on  the  geometrical 
aspects  of  the  Einstein  equations,  and  to  avoid  unnecessary 
complications  the  material  sources  have  been  considered  in  a 
simple  way.  The  dust  or  perfect  fluid  idealizations  have 
been  the  standard  models  for  matter  and  electromagnetic 
effects  have  been  treated  by  considering  charged,  unpolar¬ 
ized  fluids.  With  these  assumptions  there  is  a  clearcut 
answer  to  the  "localization"  question:  how  does  one  express 
the  total  energy-momentum  distribution  in  terms  of  material 
and  field  variables?.  One  simply  sums  the  two  standard  pres¬ 
criptions  for  material  energy -moment urn  and  field 
energy-momentum. 
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The  motivation  for  a  fully  covariant  formulation  of 
polarization,  spin,  elasticity,  etc.,  comes  from  the 
increasing  desire  to  describe  all  physical  theory  in  the 
framework  of  general  relativity.  Many  of  the  astrophysical 
objects  currently  under  observation  display  both  interesting 
physical  phenomena  (high  pressures,  large  internal  energies) 
and  intense  gravitational  fields.  For  example,  neutron  stars 
are  thought  to  possess  thick  solid  crusts  and  large  magnetic 
fields  that  are  disturbed  easily  by  the  smallest  distortions 
of  the  crust  (neutron  starquakes ) [ 1 ] .  A  knowledge  of  the 
total  energy -moment urn  of  such  a  star  forms  the  starting 
point  for  discussions  of  the  star's  gravitational  field. 

This  thesis  provides  answers  to  the  following 
questions:  how  may  spin  be  incorporated  into  gravitation 
theory?  What  form  does  the  electromagnetic  energy -momentum 
tensor  take  for  polarized  material?  It  must  be  stressed 
immediately  that  these  are  not  the  only  answers  to  the 
problems,  but  they  do  have  the  advantage  of  simplicity  over 
other  possibilities. 

Recent  work  on  the  second  question,  the  long  out¬ 
standing  " Abraham-Mi nkowski "  controversy  [2-6],  has  led  to 
an  understanding  that  the  question  is  not  really  well  posed 
and  affords  no  truly  unequivocal  solution.  The  various 
proposals  amount  to  different  allocations  of  interaction 
energy-momentum  between  matter  and  field.  Fortunately  for 
the  general  relativist,  the  fundamental  quantity  that 
couples  to  the  gravitational  field  is  the  total 
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energy -momentum  which  is  independent  of  the  particular 
convention  adopted  for  its  split  into  matter  and  field 
parts.  In  chapters  3  and  5  a  fairly  natural  splitting 
suggests  itself  and  is  adopted  for  the  purpose  of  phenomen¬ 
ological  description. 

Let  us  now  consider  our  first  question,  concerning  the 
gravitational  effects  of  spin.  Spacetime  theories  generally 
consist  of  a  geometrical  model  for  spacetime  together  with  a 
set  of  field  equations  relating  the  geometry  to  material 
sources.  In  particular,  Einstein's  general  relativity 
assumes  spacetime  to  be  a  four-dimensional  Riemannian 
manifold  with  the  Einstein  field  equations  coupling  the 
Einstein  tensor  to  a  symmetric  tensor  describing  the  total 
( non-gravi tat ional )  energy-momentum.  Because  of  the  general 
acceptance  of  Einstein's  theory,  generalizations  that  take 
spin  into  account  are  not  likely  to  be  considered  unless 
they  reduce  to  conventional  general  relativity  in  the 
absence  of  spin.  Does  Einstein's  theory  need  to  be  super¬ 
seded  by  a  more  general  theory? 

Consider  a  medium  whose  spin  is  of  orbital  (non- 
-quantum)  origin,  for  example  a  gas  of  spinning  particles 
where  each  particle  is  an  extended  body  such  as  a  rotating 
star  or  galaxy.  General  Relativity  was  originally  developed 
as  a  theory  valid  on  a  large  rather  than  atomic  scale,  in 
the  sense  that  it  is  used  for  for  continuous  matter  dist¬ 
ributions.  It  is  therefore  natural  to  use  conventional  gen¬ 
eral  relativity  for  both  the  microscopic  gravitational  field 
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(that  varies  sharply  as  one  moves  between  and  through  the 
constituent  particles)  and  the  macroscopic  gravitational 
field  generated  by  the  bulk  properties  of  the  medium.  The 
gravitational  effects  of  spin  would  be  inferred  from  the 
dependence  of  the  macroscopic  total  energy -moment urn  tensor 
on  spin. 

Some  general  remarks  may  be  made  that  indicate  how  spin 
will  contribute  to  the  macroscopic  total  energy-momentum 
tensor.  Starting  from  the  energy -momentum  distribution  of  an 
extended  body,  we  find  it  convenient  to  describe  the  evol¬ 
ution  of  the  body  in  terms  of  its  total  four -momentum,  a 
central  point,  and  the  angular  momentum,  quadrupole  and 
higher  moments  defined  with  respect  to  this  "centre".  This 
enables  us  to  visualise  the  motion  of  the  body  in  terms  of 
its  overall  motion  (motion  of  the  central  point)  and  the 
motion  relative  to  the  centre.  If  the  central  point  is 
reasonably  well  chosen  we  may  say  that  the  actual  body  and 
its  development  is,  in  some  sense,  given  to  a  first  approx¬ 
imation  by  the  centre  and  its  motion.  The  orbital  angular 
momentum  (spin)  and  the  multipole  moments,  summarizing  the 
body's  detailed  internal  structure,  provide  higher  levels  of 
approximation.  In  the  same  manner  a  distribution  of  electric 
charge  and  current  is  represented  by  the  total  charge,  total 
current,  electromagnetic  dipole  moments  and  higher  moments. 
It  seems  clear  that  one  might  decompose  any  distribution 
into  "gross"  and  "polarization"  parts  by  definition  of 
suitable  multipole  moments.  In  the  same  way  that 
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electromagnetic  dipole  moments  are  associated  with  a 
polarization  electric  four-current ,  there  should  exist 
"polarization"  energy -momentum  due  to  spin,  consisting  of 
the  divergence  of  a  third  rank  tensor  formed  out  of  the 
spin.  In  relativity  theory  the  four -momentum  and  four- 
-velocity  of  a  particle  are  not  parallel  even  for  the  most 
"appropriate"  choice  of  central  point.  This  implies  that  the 
gross  and  polarization  parts  of  the  total  symmetric  energy- 
-momentum  tensor  will  both  be  asymmetric. 

We  have  reasoned  that  since  orbital  spin  is  merely  a 
part  of  the  multipole  description  of  energy-momentum,  it 
should  influence  the  macroscopic  geometry  only  in  so  far  as 
it  determines  a  polarization  contribution  to  the  total 
energy -momentum.  Quantum-mechanical  spin,  however,  is  an 
elementary  notion  not  reducible  to  that  of  energy-momentum. 
It  may  be  argued1  that  it  should  perhaps  have  a  geometrical 
significance  of  its  own,  both  (canonical)  energy-momentum 
and  spin  flux  independently  leaving  an  "imprint"  on  the 
spacetime  geometry.  This  is  a  feature  of  the  "torsion" 
theories  of  Weyl  [33],  Sciama  [34] 2  and  Kibble  [35],  dev¬ 
eloped  further  by  Trautman  [36],  Hehl  [37]  and  others  [38]. 

Torsion  theory  has  developed  mainly  out  of  local  gauge 
theory  for  the  Poincare  group  and  not  from  any  incompat¬ 
ibility  of  conventional  general  relativity  and  spin.  One  may 


1 Cf .  [40]  and  [38,  page  394].  The  latter  reference  contains 
a  general  overview  of  torsion  theory  and  an  extensive 
reference  list. 

2Sciama  has  also  suggested  the  use  of  an  asymmetric  metric 
tensor  in  connection  with  spin  [39]. 
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just  as  easily  assume,  therefore,  that  spin  of  quantum 
origin  and  spin  of  orbital  origin  both  influence  the  space- 
time  geometry  in  the  same  way.  One  would  then  retain  the 
standard  form  of  Einstein's  theory  with  the  total  spin  of 
both  types  determining  the  polarization  energy -moment urn. 

When  the  field  equations  of  torsion  theory  are  translated 
into  an  equivalent  Riemannian  form  they  differ  from  those  of 
conventional  general  relativity  by  a  term  quadratic  in  the 
spin  on  the  right-hand  side  of  the  usual  Einstein  field 
equations.  At  the  ordinary  macroscopic  level  where  the  grav¬ 
itational  effects  of  spin  are  insignificant  [41]  both  the 
Riemannian  theories  and  the  non-Riemanni an  theories  lead  to 
the  usual  (spin-free)  form  of  the  energy -moment urn  tensor. 

The  dynamics  of  polarized  media  will  be  formulated  in 
this  thesis  entirely  within  the  framework  of  standard  gener¬ 
al  relativity.  For  orbital  spin  this  is  the  obvious  choice, 
for  elementary  spin  it  has  the  advantage  of  simplicity  over 
generalizations  of  Einstein's  theory. 

In  chapter  two  the  general  form  of  the  Lagrangian 
equations  of  motion  for  a  spinning  multipole  particle  in 
given  external  fields  is  derived  from  an  action  principle 
involving  an  unspecif ied  Lagrangian.  The  special  case  of 
motion  in  external  electromagnetic  and  gravitational  fields 
is  considered. 

Chapter  three  concerns  the  Lagrangian  dynamics  of 
continuous,  polarized  media  with  intrinsic  spin.  The  Ein¬ 
stein  field  equations  are  derived  together  with  the 
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non-gravi tat ional  field  equations  and  the  material  balance 
laws  for  momentum  and  spin.  A  simple  method  of  splitting  the 
total  energy-momentum  tensor  into  matter  and  field  parts  is 
proposed.  For  an  electromagnetic  field  this  leads  to  a 
"generalized  Minkowski"  electromagnetic  energy-momentum 
tensor . 

In  chapter  three  the  dependence  of  the  Lagrangian 
density  on  derivatives  of  the  fields  is  restricted  to  first 
derivatives  in  order  to  simplify  the  discussion.  Chapter 
four  considers  the  general  case  where  the  Lagrangian  density 
includes  higher  derivatives  of  the  fields.  The  partition 
into  matter  and  field  produces  an  electromagnetic  tensor 
depending  on  all  derivatives  of  the  electromagnetic  field 
and  all  multipole  moments. 

Chapter  five  examines  the  concept  of  gravitational  pol¬ 
arization  (polarization  of  energy-momentum)  by  comparison 
with  the  traditional  multipole  analysis  of  electromagnetic 
polarization.  Using  the  multipole  formalism,  the  particle 
equations  of  motion,  balance  laws  and  general  structure  of 
the  total  energy -momentum  are  rederived.  Quantities  such  as 
four-momentum  and  spin,  whose  interpretation  in  the  Lagrang¬ 
ian  theory  is  unclear,  acquire  a  direct  physical  meaning. 

In  chapter  six  we  use  a  Lagrangian  technique  of  calcul¬ 
ation  in  order  to  develop  the  multipole  formalism  in  curved 
spacetime.  We  consider  the  expansion  of  the  Lagrangian  of  a 
simple  (non-spinning)  particle  and  then  show  that  variation 
of  the  particle  world-line  results  in  the  spinning  particle 
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equations  when  the  expanded  form  of  the  Lagrangian  is  used 
in  the  action  principle.  Expressions  for  the  four -momentum 
and  spin  of  an  electromagnetic  multipole  in  a  Einstein- 
-Maxwell  field  are  obtained.  They  are  shown  to  be  the  same 
as  those  proposed  by  Dixon  [11]. 


- 


2.  SINGLE  PARTICLE  EQUATIONS  OF  MOTION 


2.1  Discussion 

Spinning  particle  equations  of  motion  have  tradition¬ 
ally  been  derived  either  from  an  action  principle  or  the 
multipole  formalism  for  an  extended  body.  The  multipole 
formalism  [9,10,11,5]  involves  integration  of  the  conser¬ 
vation  law  for  the  total  energy  tensor  over  space! ike 
sections  of  the  particle  world  tube1.  With  a  given  material 


energy  tensor  T 


a  6 


one  must  choose  a  set  of  tube  sections 


(mat) 

E ,  a  central  world  line  Za(x)  and  some  method  of  propagating 
the  material  four-momentum  density  T°^  .  n  2  from  each  point 

of  E  to  the  point  of  intersection  of  Za(r)  with  e  (reducing 
to  parallel  propagation  for  zero  curvature) .  In  special 
relativity,  any  given  choice  of  E' s  determines 


a 


P  00  = 


-  J 

E 


T 


aS  n 

(mat)  3 


dE 


(2.1  ) 


Any  choice  of  Z“(t)  parametrizes  the  z' s  and  defines3 


S“B(z(x)  ,  Z  (x  ) ) 


2  /  (x[ol-Z[a)  TB]y  ,  n  dz  .  (2.2) 

*,  (mat )  y 


Different  choices  of  central  world- line  and  tube  sections 
lead  to  different  so-called  "auxiliary  conditions",  such  as 
Sa^uD  =  0,  Sa^p0  =  0  and  others  [13].  In  the  general 

p  p 


^or  a  discussion  see  Ehlers  [12]. 

2Where  ng  is  the  unit  timelike  normal  vector  field  of  E. 

3  The  four-momentum  and  spin  of  the  particle  are  pa  and  Sa& 
plus  terms  depending  on  the  applied  fields. 
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relativistic  generalization  of  (2.1)  and  (2.2)  the  most 
convincing  arguments  for  the  choice  of  propagation  law  for 
I  (mat)  n3  have  been  given  by  Dixon  [11].  Interestingly,  the 
bitensor  field  he  uses  to  carry  .  n  to  Za  in  his  defin- 
ition  of  spin  differs  from  the  one  used  in  his  definition  of 
fourmomentum. 

Obviously  the  dynamical  quantities  momentum,  spin,  etc. 
depend  upon  the  specific  choice  made  for  space  sections, 
central  world-line  and  propagation  law.  When  a  choice  is 
made  at  the  outset  of  the  derivations  it  becomes  difficult 
to  answer  the  question:  to  what  extent  is  the  form  of  the 
equations  of  motion  independent  of  the  choice?  It  is  easier 
for  Lagrangian  dynamics  to  provide  an  answer  than  the 
multipole  formalism  because  it  is  fairly  simple  to  obtain 
the  equations  of  motion  from  a  variational  principle.  It 
will  become  clear  in  the  next  section  (and  also  in  chapters 
5  and  6  from  a  different  viewpoint)  that  the  equations  are 
very  general  in  nature.  This  generality  may  be  seen  in 
Dixon's  equations  (7.1),  (7.2),  [11c].  His  careful  analysis 
does  not  rely  on  a  specific  choice  of  Za(x).  For  tube 
sections  Dixon  uses  "geodesic"  hypersurfaces  orthogonal  to 
an  arbitrary  timelike  vector  field  along  Za.  Dixon  only 
introduces  specific  (and  natural)  conventions  for  Za  and  £ 
for  convenience  at  a  later  stage. 

We  derive  the  general  form  of  the  Lagrangian  equations 
of  motion  for  a  spinning  particle  having  arbitrary  multipole 
structure  in  arbitrary  external  fields  (eqs.  (2.22)  and 


(2.23))  and  then  examine  the  form  of  the  equations  for 
motion  in  a  Maxwel 1 -E i nstei n  field.  The  treatment  is  inde¬ 
pendent  of  special  interaction  models  in  that  we  do  not 
specify  the  functional  form  of  the  Lagrangian,  but  only  the 
variables  on  which  the  Lagrangian  depends.  This  generality 
of  the  Lagrangian  also  makes  a  priori  commitment  to  part¬ 
icular  conventions  such  as  a  definite  auxiliary  condition 
unnecessary.  The  resulting  equations  represent  the  canonical 
substructure  to  which  every  detailed  model  of  particle 
motion  must  conform.  The  two-page  calculation  below  is  a 
concise  derivation  of  spinning  particle  equations  avoiding 
the  complicated  workings  of  the  multipole  formalism.  (Of 
course,  there  is  little  physical  insight  in  definitions  such 
as  Pa  =  3  L / 3  va  compared  with  (2.1).)  Section  2.4  reviews  the 
special  case  of  a  multipole  particle  in  external  electromag¬ 
netic  and  gravitational  fields. 

Action  principles  have  been  used  previously  by  Frenkel 
[14]  with  choice  of  a  special  Lagrangian  and  more  generally 
by  Barut  [15,  page  77]  without  specifying  the  functional 
form  of  the  Lagrangian.  They  have  derived  special  relativ¬ 
istic  equations  of  motion  for  spinning  particles  in  electro- 
magnetic  fields,  both  assuming  S  ug  =  0 .  In  general  rela¬ 
tivity,  Kunzle  and  others  [16,17]  have  obtained  equations 
for  spinning  dipoles  in  Einstein-Maxwel 1  fields  from  a 
Lorentz- i nvar i ant  pre-symplectic  1-form:  they  adopt  the 
auxiliary  condition  S  p  =  0.  Fuchs  [18]  has  used  an 

p 

approach  similar  to  ours  with  tetrad  variables.  Yet  other 


- 


■ 


derivations  of  particle  equations  of  motion  proceed  on  a 
more-or-less  ad  hoc  basis  [19]. 
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2.2  Definitions  and  General  Identities 

We  begin  by  introducing  some  general  formulae  for 
relative  tensor  fields  $  (xa),  and  then  develop  the  identit- 
ies  (cf.  Belinfante  and  Rosenfeld  [7])  which  flow  from  the 

condition  that  a  function  L($  )  be  a  (relative)  scalar1. 

A 

A  set  of  relative  tensor  fields  will  be  denoted  in  a 
compact  way  as  $  ,  i.e.,  capitalized  Latin  indices  denote  a 
set  of  tensor  labels  and  the  respective  tensor  component 
indices.  Thus, 


(<f>  ) 

a  a ....  a 
1  m 


a  .....  a 
m+1  n 


a  .  =  1  to  4 
i 

3.  1 , 2  j  .  .  . 


where  m  and  n  depend  on  the  label  a.  Let  $A  denote  the  same 
set  of  tensors  with  covariant  component  indices  raised  and 
contravar i ant  component  indices  lowered: 


$ 


A 


a 


(4>  ) 

a 


1  *  ' 


a 

m 


a 


m+1 


Sometimes  an  index  A  is  underlined  with  a  tilde  ~  as  a  means 
of  indicating  that  A  and  A  have  different  1 abel -ranges .  For 

example  y  below  is  used  to  denote  the  set  of  tensors  ip .  and 

A  A 

their  first  covariant  derivatives:  y .  e  ( ^  A ,  ip  .  ).  A 
repeated  index  implies  summation  over  the  respective  tensor 


^ore  details  may  be  found  in  [8]. 


I 
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indices  and  over  the  label -range  concerned. 

Under  the  co-ordinate  transformation  =  xa(x^)  the 
relative  tensor  $  transforms  linearly: 

%(x)  =  AaB(xPcj)*b(x)  *  xP0  s  3xp/axa  .  (2.3) 

g 

(Of  course,  AA  vanishes  unless  a  =  b  so  the  summation  over 
B  in  (2.3)  does  not  couple  different  tensors.)  The  infinit¬ 
esimal  generators  of  the  transformations  (2.3)  are 


(i 


CT 


P 


( 3  A  B/ 3  XP 
A  a 


(2.4) 


(The  exact  form  of  these  is  given  in  Appendix  1.)  Explicit 
construction  of  these  generators  follows  easily  from  the 
recursion  formulas 


(I  B3)  a  =  (I  B)  a63  -  6B6a63 

Aa  p  A  p  a  Aap 


(2.5) 


,  T  a  B  N  a  B .  a  a  Baa 

(IAB)P  =  (IA  A  5B  +  SASPSB 


(2.6) 


and  (I)  a  =  -  w6°  (2.7) 

p  p 

for  a  relative  scalar  <(>  of  weight  w  (<j>(x)  =  <j)  ( x )  |  ax/ ax  |  w)  . 

T) 

Here  6A  =  1  if  a  =  b  and  the  respective  tensor  indices  are 

•D 

equal,  and  6A  =  0  otherwise. 

Useful  formulas  which  will  be  required  in  the  sequel 


are 


■ 
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$ 

A 


T 


d  $.  +  T9  (I  B)  % 
t  A  ax  A  p  B 


(2.8) 


$  i  —  $  i 

A  |  y  v  A  |  vy 


Rp 

ay  v 


d.B) 


(2.9) 


where  the  stroke  denotes  covariant  differentiation  with 
respect  to  a  symmetric  affine  connexion.  From  (2.3)  and 
(2.4)  the  change  of  $A  under  the  infinitesimal  co-ordinate 
transformation  xp  =  x9+e p ( x )  is 


$  (x)  -  $.(x)  =  (I  B)  a$(3  Ep)  .  (2.10) 

A  A  A  p  rS  a 

From  (2.10)  we  derive  at  once 


aL  (I,?)  V  +  wS°L  =0  (2.11) 

-  A  p  B  p 

B'F 

A 

as  the  condition  that  the  function  L(v  )  of  the  relative 

A 

tensors  'F  transforms  as  a  relative  scalar  of  weight  w. 

A  more  restrictive  identity  can  be  derived  from  (2.11) 

in  the  case  of  a  scalar  density  (w=1)  l_U  ,  $  ,  )  depending 

A  A  |  p 

on  a  set  of  relative  tensor  fields  and  their  first1 

A 

covariant  derivatives.  We  define  the  variational  derivative 


SL/Sip,  =  LA  -  LAa,  (2.12) 

A  |  Ot 

where 

LA  h  dL/dipA  ,  LAa  =  8L/^A|a  (2.13) 


1  The  case  where  higher-order  derivatives  occur  is  not 
needed  immediately  and  its  consideration  is  deferred  to 
chap . 4  f  sec . 1 . 
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and  the  tensor  densities 


V  5  LSP  -  *a|plA°  •  UTOp  H  LAt(IaB)p%b  .  (2.14) 


With  ^  A  =  ( ^  A ,  tPAi)p  condition  (2.11)  is 


lA  Op°*B  +  «-Aa(IAoBB)p0*B|S  +  ^  ° 


(2.15) 


Further,  we  note  from  (2.5)  that 


LA“(t  BB)  %  .  =  LAo(l.B)  ,  -  LAV, 

Aa  pB|3  ApB|a  A  p 


( LAa  (  i  B )  %  ),  -  LAa,  (t  B)  <V  -  LAV, 

A  p  B  |  a  a  A  p  B  A  p 


which  inserted  into  (2.15)  gives 


U 


Ta  ,  +  t  0  +  sL  (I  B)  =  0 

p  |  x  p  —  A  p  B 


(2.16) 


6  ip 


valid  for  an  arbitrary  scalar  density  L  ( ip  ,  \p  i  ). 

A.  n  01 

In  equations  (2.12  to  16)  the  letter  L  printed  bold 


face  denotes  a  scalar  density.  Tensor  densities  t  a,  UTa 

p  p 

and  SpTa  formed  from  L  are  also  printed  bold  face. 


2.3  Spinning  Multipole  in  Given  External  Fields 

Let  xy  =  xy(t)  be  the  equation  of  the  particle  world¬ 
line  in  terms  of  an  arbitrary  scalar  parameter  t,  x  the 
proper  time,  and 


16 


uy  =  dxy/dr  ,  vy  =  dxy/dt 


the  normalized  and  unnormalized  four -ve loci t ies .  The  spin  of 
the  particle  is  described  by  the  gyration  of  an  orthonormal 
tetrad  e  (t)  defined  on  the  world  line: 

a 


(a)  (b) 

n  ,ev  ;e:  =  g  Q 

aba  p  a p 


n  u  =  nab  =  diag(l  (a,b  =  1  to  4) 

a  b 


(a)  (b)  =  ab 

a  3 


2.17) 


The  equations  of  motion  are  assumed  to  be  derivable 
from  a  par  ameter- i  nvar  i  ant  action  principle  s/Ldt  =  0,  for 

variations  of  xu ( t )  and  e^a^(t)  with  fixed  endpoints  xy ( t .  )  , 

a  l 

(a  ) 

e  (t.)  (i  =  1,  2).  (The  parameter  invariance  eliminates 
the  constraint  u^uU  =  -1,  see  [8,  p . 1 1 ] . )  The  Lagrangian  is 
an  unspecified  scalar  function 

L  =  L(vy,  e(a),  4(a),  »  )  (2.18) 

a  a  A 

where  e^  =  Se^Vst  (absolute  derivative)  and  the  set  of 

a  a 

tensors  $  comprises  the  external  fields  <f>  ,  the  Riemann 

A 

tensor  RV  .  and  their  covariant  derivatives.  We  have 

3  y<5 

assumed  that  explicit  dependence  of  L  on  gag  has  been 
eliminated  by  the  use  of  (2.17). 

We  define  the  canonical  momentum  ,  the  spin  angular 


' 
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momentum  S 


pa 


=  S 


[pa] 


and  the  multipole  moments  M~  by 


P 

y 


2 


e 


(a) 

[p 


3  L 


3  G 


TsT) 

a] 


(2.19) 


(The  derivatives  3L/3$A  =  are  defined  such  that  dL  = 

A  A 

M~d$A  and  M~  has  the  same  algebraic  symmetries  as  $  .  For 
example,  if  L  =  F^Xyv  with  F  ^  antisymmetric,  then  L  = 

F  XC^]  and  3  L / 3  F  is  defined  as  X^-yv^  .) 

y  v  ]i  v 

The  equations  of  motion  for  the  spin  Sp°  are  obtained 
by  variation  of  e^a^(t)  with  fixed  end-values,  holding  the 
world-line  fixed.  Since  in  this  section  we  treat  the  space- 
time  geometry  as  prescribed,  equations  (2.17)  impose  10 
constraints  on  the  16  variations  6e^  .  The  resulting 

equations  of  motion  reduce,  after  eliminating  the  Lagrange 
multipliers  [8,  page  10],  to  the  six  equations 


5  L 


(a)  _ 


6  e 


(a  )  p  ] 
[a 


=  0 


6  L 


3L 


6  e 


(a) 


a 


3  e 


(a) 


a 


&  t 


3  L 


3  e 


(a) 


(2.20) 


On  the  other  hand,  the  identity  (2.11)  with  w  =  0  gives 


3  VP 


3  L 


g  (a )  +—(I?}  a$ 

3  g(a )  P  9V  ^  p  5 

a 


0  . 


(2.21 ) 


Differentiating  the  second  of  (2.19),  we  find  with  the  aid 
of  (2.21)  that  eq.  (2.20)  is  equivalent  to 


^he  parameter  invariance  of  I  implies  that  Py  and  Sp  are 
both  independent  of  parametrization.  is  proportional  to 
dx/dt  and  can  be  fixed  by  setting  t  =  t  in  the  equations  of 
motion.  The  definition  of  S  a  is  discussed  in  appendix  10 
and  in  [8,  p.67].  p 


- 


. 
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-  6Spa/6t  =  P[pva]  +  M~  (I  ?)  [pa-*  $ 

2  a  b 


(2.22 


The  equations  of  motion  for  the  linear  momentum  are 
obtained  from  an  infinitesimal  displacement  of  the  world- 
-line,  holding  e^  J  fixed  by  parallel  propagation.  We 
consider  a  1 -parameter  family  of  time- like  curves  xu(t,e) 
with  orthonormal  tetrads  e_^a^(t,e)  defined  on  them,  and 

extremize  1(e)  =  -J4‘  ■c- - 1 

(i  =  1,  2),  Se^Vse  =  0.  This  yields,  for  arbitrary 

a 

var i at  ions  8XU/3e , 


2|_dt  subject  to  fixed  xy(t.,e)  =  xy(t.,0) 

1 1  i  1 


3  L  6  Vy  3  L 
3  vy  6  e 


3  e 


(a) 


6  £ 


6  e 


(a) 


t  5t 


,  3L  , 
+  —  $ 


3  X' 


3$  A.|  u  3e 

A 


dt  =  0 


Noting 


u 

'  3  xy ' 

6 

fse  (a)l 

6 

6 

a 

6  £ 

6  t 

[3£  j 

6  £ 

l  6 1  . 

=  e(a)Rp  ,  vX3xu 

p  o\u  ^ 


and  integrating  by  parts  gives  the  equations  of  motion 


6P  /fit  =  . 

U  2  p  a  Xy 


SpavX  + 


A  ha 


(2.23) 


2.4  Motion  in  a  Maxwel 1 -Einstein  Field 

It  will  be  convenient  to  use  the  following  notation:  a 
tilde-underlined  index  a(n)  denotes  the  symmetrized  set  of 
indices  {a  a(0)  denotes  the  empty  set.  Repeated 

indices  imply  summation.  Thus, 


^Aa  (n  )  y 


Ba  (n  ) 


.X 


xay 

A  ( ot j  •  •  a  )  y 


B 


B  (a 


i  f  n  =  0 


s  if  n  =  1,2, ..J  (2.24) 

X  •  • an  > 
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As  a  special  case  of  Section  2.3,  let  the  set  of 
tensors  $A  comprise  the  four-vector  potential  A  ,  the  elect¬ 


ee 


romagnetic  field  tensor  F  „  =  23 ,  Ao1,  the  Riemann  tensor 

a  3  la  B  J 


ct  * 

R  .  .  and  the  symmetrized 1  covariant  derivatives  of  F  „  and 

3y6  aB 


R 


a 


3  Y  6 


,a 


* ,  =  (A  ,  F.  ,  .  ,  R‘.  n 


a 


By  j  A  (n)  ’  By 6  |  A  (n)  * 


=  0,1,2 , ... )  .  (2.25) 


Define  the  electromagnetic  multipole  moments 


s  2  8L/3F  |x(n)  (t-T) 


(2.26) 


.n  +  2 


and  the  2  -pole  gravitational  moments1 2 


By  5  e 


-n 


a 


(x)  =  2  3L/3R 


a 


3 y<5  |  e  (n) 


(  t  =  T  ) 


(2.27) 


We  assume  the  (undifferentiated)  potential  Aa  enters  L  only 
through  a  bilinear  interaction  term  eA^v01  and  define  a  gauge 
invariant  "Kinetic  momentum" 


1  This  is  not  a  restriction,  since  an  unsymmetrized  covar¬ 
iant  derivative  can  always  be  reduced  to  a  set  of  the  form 
(2.25)  with  the  aid  of  the  Ricci  commutation  relations 
(2.9) . 

2  This  corresponds  (apart  from  a  numerical  factor)  to  the 
"reduced"  moment  integral  Y<5  e  1  *  *  *  en  defined  for  extended 
mass  distributions  by  Dixon  [lie,  eq.  (5.33)].  (2.26) 
corresponds  to  [lie,  eq.  (5.32)].  With  the  convention  that 
Ma  =  3L/3$a  share  the  same  symmetries  as  $A,  (2.26)  is 
antisymmetric  in  B,  y,  and  symmetric  in  its  last  n  indices. 
(2.27)  has  the  algebraic  symmetries  of  the  Riemann  tensor  in 
its  first  four  indices  and  is  symmetric  in  its  last  n 

i ndices . 
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p  =  P  -  eA 

a  a  a 


(2.28) 


With  (2.25  to  2.28),  the  translational  equations  of  motion 
(2.23)  reduce  to 


6p  /6t  =  -  \  R  .  _  uBSy<S  +  eF  „uB 
a  Z  a  3  y  o  a  3 


(2.29) 


i  ” 

+  7  l  F( 


2 Q1  3  Y  |  X  (n)  a 


m 


3y  X  (n) 


v  3 

J0R  y6X  |  e  (n  )  a  ^  3 


y 5  X  e  (n) 


The  spin  equations  of  motion 1  reduce  to 


\  6SaB/6T  =  p[aue] 


(2.30) 


V  c  [« 

n=0  Yl^(n) 


mB]yX(n)  _  1  £  +  mYVX(n)B] 

2  _n  yy  (X  (n)  •  ) 

n=  U  ~ 


oo 

4  l  R[“ 


n=  0 


ou 

q8]yuvX(n)_  (n+1)R« 
n=  0 


yiiV  |  X  (n)  V|  ~  YU v  i  (X(n)»)^6 


yy  vX  (n) 3  ^ 


(  A  dot  placed  above  an  index,  together  with  a  square  bracket,  denotes 

[ .  .  ] 

antisymmetrization  of  that  index.  Thus  AaB(3  =  A^aB^  .) 


’See  Appendix  1. 
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3.  GENERAL  RELATIVISTIC  FIELD  EQUATIONS  FOR  POLARIZED  MEDIA 

WITH  INTERNAL  SPIN 


3.1  D i scuss i on 

In  this  chapter  we  derive  the  Lagrangian  dynamics  of  a 
polarized,  continuous  medium  with  internal  spin,  in  inter¬ 
action  with  fields  generated  by  the  medium. 

Reviewing  the  existing  literature,  we  find  the  main 
approach  in  describing  the  interaction  of  elastic  materials 
or  fluids  with  a  Maxwell  or  E i nstei n-Maxwe 1 1  field  has  been 
a  variational  one  [20,  21,  22,  4,  23,  24].  Of  these  refer¬ 
ences  ,  in  a  general  relativistic  framework,  Schopf  has 
considered  dielectric  fluids  in  [22]  and  Maugin  discusses 
magnetized  elastic  media  in  [4a]  and  magnetohydrodynamics  in 
[4b].  Maugin  and  Eringen  have  used  Lorentz  invariance  to 
derive  balance  laws  for  for  polarized  elastic  solids  with 
electronic  spin  in  special  relativity  [23].  Other  references 
and  a  general  overview  of  the  Lagrangian  dynamics  of  contin¬ 
uous  media  can  be  found  in  Soper's  book  [24]. 

Carter  and  Quintana  [25]  have  considered  the  general 
mathematical  formalism  for  continuous  media  in  general  rela¬ 
tivity.  They  discuss  the  covariant  formulation  of  elasticity 
theory,  the  concept  of  perfect  elasticity,  and  they  mention 
the  possible  applications  to  the  areas  of  gravitational  wave 
detection  and  astrophysics  ( Ruderman  [1]  has  emphasized  the 
probable  crystalline  nature  of  the  outer  regions  of  neutron 
stars )  . 
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Fluids  with  internal  spin  have  been  discussed  by 
several  authors  [26].  It  is  well  Known  that  the  presence  of 
spin  leads  to  the  asymmetry  of  the  material  energy-momentum 
tensor.  (For  polarized  matter  the  material  energy -moment urn 
tensor  is  generally  asymmetric  even  without  spin,  when 
certain  matter-field  interaction  momenta  are  allocated  to 
the  mater i a  1 . ) 

The  present  analysis  is  based  on  a  Lagrangian  density 
which  is  an  unspecified  function  of  three  scalar  fields  am 
(m  =  1,2,3)  (co-moving  parameters  of  the  medium),  an  ortho- 

normal  tetrad  field  e  ( represent i ng  a  set  of  spin  axes) 

a 

and  arbitrary  field  variables.  Physically,  this  represents  a 
general  continuous  medium  with  spin  interacting  with  the 
external  fields  treated  as  self-consistent  backgrounds.  The 

tetrad  also  serves  to  specify  the  metric  field.  Variation 

C  3-  ^ 

with  respect  to  the  sixteen  tetrad  components  ea  yields 
six  equations  for  the  spin  angular  momentum  and  the  ten 
gravitational  field  equations.  The  translational  equations 
of  motion  of  the  material  are  then  obtainable  (modulo  the 
non-gravi tat iona 1  field  equations)  either  from  the  con¬ 
tracted  Bianchi  identities  or  from  "variation  of  the  world 
lines".  Stress  contributions  to  the  material  energy-momentum 
tensor  arise  from  the  dependence  of  the  Lagrangian  upon  the 
gradients  3aam  (cf.  [24,  page  57]). 

In  section  6  the  analysis  is  applied  specifically  to  a 
medium  interacting  with  an  Einstein-Maxwel 1  field  to  give 
clear  formal  expression  to  the  general  relativistic 


' 
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extension  of  Lorentz  dielectric  theory.  The  specialization 
to  non-viscous1  spinning  fluids  and  "dust"  is  discussed  in 
section  7. 

3.2  Kinematics  of  the  Continuum 

Before  deriving  the  system  of  field  equations  we 
discuss  the  general  Kinematics  of  the  medium.  In  terms  of  a 
"numerical"  flux  vector  field  Ny(x)  (particle  flux)  the 
number  density  of  particles  per  metric  volume  n(xy),  meas¬ 
ured  in  a  local  rest  frame  at  xy ,  is  defined  by 

Ny  =  nuy  ,  uyu  =  -1  .  (3.1) 

y 

The  motion  of  the  medium  is  described  by  the  world  lines, 
defined  by  the  equations 

3xy(am,t)/3t  =  vy  =  uyd-r/dt  (3.2) 

as  integral  curves  of  Ny.  The  three  am  are  "co-moving" 
parameters  (uy3^am  =  0)  and  t(t)  is  an  arbitrary  parameter 
along  the  curves,  t  being  the  proper  time  along  each  curve. 
For  an  elastic  solid  one  may  think  of  the  co-moving  param¬ 
eters  am  as  "attached"  to  each  lattice  point,  each  lattice 
point  having  associated  with  it  fixed  values  of  a™  through¬ 
out  its  motion.  A  fluid  would  have  its  macroscopic  motion 


'The  non-di ss i pat i ve  nature  of  the  dynamics  is  an  assumption 
inherent  in  the  Lagrangian  treatment. 
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described  by  "stream  lines",  each  labelled  by  a  particular 
set  of  values  for  am.  For  both  fluids  and  solids  the  number 
of  particles  in  an  infinitesimal  flux  tube  d3a  is  assumed  to 
be  a  constant  of  the  motion  given  by  N(am)d3a.  In  terms  of 

the  gradients  3  am  and  N(am)  the  numerical  flux  may  be 

y 

written  as 

Nu(x)  =  N(am)nya6Y(3  a1)  (3  a2 )  (3  a3)  (3.3) 

0*  P  y 

where  is  the  completely  antisymmetric  permutation 

tensor1.  Conservation  of  particle  number  implies  that 


Defining  v  by 

mn 


ran 


Y 


mn 

Y 


gy  v  ( 3  am)  ( 3  a11) 

3  y  '  '  v 


(3.5) 


distances  between  neighbouring  points  am  and  am+dam  of  the 
material  in  the  local  rest  frame  are  given  by 


(ds2)^  =  ymndamdan  .  (3.6) 

A  useful  form  for  the  number  density  n  is  its  expression  in 
terms  of  N(am)  and  Ymn, 


i.eyoteY  =  /^g’nya^y  where  eyaBy  are  the  Levi-Civita 
permutation  symbols. 
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n(xy)  =  N ( am) (de tYmn) ^  .  (3.7) 

The  material  will  be  assumed  to  be  in  adiabatic  motion  with 
the  entropy  per  particle  S  constant  along  each  world  line,  S 
=  S(ara).  The  entropy  current  Sy  =  S(am)Ny  according  to  (3.4) 
then  satisfies  Su i  =  0. 


3.3  Gravitational  Field  Equations  Derived  from  Variation  of 
Tetrad 

The  gravitational  field  equations  are  derived  from 

(a) 

variation  of  an  orthonormal  tetrad  field  ea  (x)  that 
satisfies 


(a)  e(b)  =  g 
aba  B  ya  B 


g a  B  e( a)  e(b) 

3  a  3 


=  n 


ab 


(3.8) 


n  ,  =  nab  -  diagd.1,1,-1)  (a,b  =  1  to  4) 

a  b 


(3.9) 


The  sixteen  component  tetrad  field  plays  a  dual  role:  the 
ten  symmetrized  products  nabe^')  e^  =  gag  define  the  metric 
and  determine  the  gravitational  field,  while  the  dependence 
of  the  Lagrangian  density  upon  the  six  angular  velocities 
wab  _  e(a)°t  determines  the  internal  spin  of  the 

medium.  Because  of  this  second  role  the  Lagrangian  density 
will  depend  explicitly  on  the  tetrad  field  in  addition  to 
its  dependence  via  the  metric  tensor  ga^-  One  could  have 
considered  a  Lagrangian  density  depending  upon  both  gag  and 
e  (and  then  varied  g„D  to  obtain  the  ten  gravitational 


a 


. 
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equations  and  varied  six  independent  components  of  the 
tetrad  field  to  obtain  spin  equations).  However,  since  the 

or  tnonorma  1  i  ty  conditions  g  =  n  e^  e^  are  ten  rela- 

a  3  ab  a  3 

tions  between  the  twenty-six  variables  (g  ,  e ^  ),  we  have 

a  3  a 

sixteen  independent  variables  and  it  is  much  simpler  to 

choose  the  e^  as  the  sixteen.  We  take,  therefore,  arbit¬ 
er 

rary  independent  variations  in  the  e^  to  obtain  ten  grav- 

a 

itational  equations  and  six  spin  equations. 

Helpful  guidance  for  the  following  calculations  comes 
from  Rosenfeld  and  Belinfante  [7],  They  consider  spin  flux 
due  to  fields  only,  while  a  total  spin  flux  due  to  both 
material  and  field  is  considered  here. 

The  spin  and  translational  equations  of  motion  for  the 
medium  and  the  system  of  field  equations  for  the  applied 
fields  are  all  obtained  from  the  four-dimensional  action 
integral 


I  = 


L(v  W 


.4 

d  x  + 


( 16  TT ) 


-1 


^  R 


d4x 


(3.10) 


in  which  we  have  made  the  specific  choice  of  the  curvature 
scalar  R  for  the  gravitational  free-field  Lagrangian1.  The 
Lagrangian  density  L  is  taken  to  be  an  unspecified  function 
of  the  fields 


According  to  Lovelock's  theorem  [27a],  any  choice 
L ( gag  ,  3eg  3 , . . . ,  3 e : • ■ • 3 e ng a 3 * • • • )  which  leads  to  second 
order  field  equations  for  minimally  coupled  sources  isa^ 
variational ly  equivalent  to  the  curvature  scalar  R  =  g  R  R 

(R  R  =  R  )• 

a3  ay3 


' 
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$ 


A 


4> 


R 


a 


gy  6 


(3.11) 


and  first  covariant  derivatives  ^  i  .  4>  A  is  an  arbitrary  set 
of  fields  interacting  with  the  gravitational  field  and  the 
medium.  We  are  restricting  the  dependence  of  L  to  only  first 
derivatives  in  order  to  simplify  the  discussion.  To  simplify 
further  we  will  also  assume  that  no  derivatives  of  R%y5 
enter.  This  is  still  sufficiently  general  to  cover  most 
cases  of  practical  interest.  (For  the  analysis  of  the 
general  case  see  Chapter  4.)  The  Lagrangian  density  L  is 

assumed  to  be  constructed  from  am  and  am,  =  3  am  so  as  to 

I  y  y 

be  independent  of  the  particular  choice  of  "material  co- 
-ordinates"  am.  (Nu(x)  and  n(x),  given  by  equations  (3.3) 
and  (3.7),  are  both  invariant  under  am  am  =  fm(an).) 

Under  an  arbitrary  variation  <se  ^  of  the  tetrad  field, 

a 

and  the  resulting  variation 


sg 


P  c 


2n  e(b)  6e(a) 

ab  ( p  a  ) 


(3.12) 


of  the  metric,  variation  6.  .L  has  contributions  from  three 

(e) 

sources : 

1.  the  explicit  dependence  of  L  on  e^  ,  eva|  ^  ; 

2.  the  variation  of  the  affine  connexion  hidden  in  the 

covariant  derivatives  e ,  <j>  .  .  ; 

a  |  3  Y  A  |  a 

3.  the  variation  of  Ra„  „  . 

gy<5 

Accordingly,  from  (2.8)  we  have  5.  ,L  = 

a  (e) 


. 
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UL/se(a))se(a)  +  UTa  srp  +  0  By66R“  .  + 

a  a  p  ax  a  pyo 


where  6L/5ij>A  and  U  were  defined  in  (2.12) 

A  p 


0  =  /=i  Q  ey6  =  3L/3R 

a  a 


a 


gy6 


and  (div)  represents  a  divergence  8  (...). 


a 


To  re-express  the  last  two  terms  of  (3. 


(a) 


$e  ,  we  note  that 


,a 


5fT„  „  =  2(sr“  ,  srp^  =  ( 6 g ) 


3  Y  6 


3  [  5  7  |  y  ] 


P 
a  x 


P  1 

(air)  2 


enable  us  to  write  the  identities 


Q  5  6  Ra  .  =  -2Q  a[yT|  STP  +  (div) 
a  3  Y  <5  P  |  y  0  r 

UTa  <srp  =  \(Usoxp  +  Spxa)  -  uT  (pa)  )  |  6 g  a 

p  a  x  2  2  l^pa 

valid  for  arbitrary  tensor  densities  Qa^y<S, 
defined  the  "spin  flux" 


SPTp 


H  2  U 


p  [  ax  ] 


Using  the  variation  in  R 

(16ir)'16(/zgR)  -  '  (  8tt  )  ”  1  6e(a)aBeg 


(div)  (3.13) 

and  (2.14), 

(3.14) 

13)  in  terms  of 

-(fig)  |p  (3.15) 

OT 

(3.16) 

+  (div)  (3.17) 

UT0  .  We  have 

p 

(3.18) 

a) 


(3.19) 


- 
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where  GaB  is  the  Einstein  tensor,  GaB  =  RaB  -  -|gaBR,  the 
action  principle  6^1  =  0  gives,  with  the  aid  of  (3.13), 
(3.16),  (3.17)  and  ( 3  .  12) 


(SnrVr^r  =  UL/«e<a))e(a)p 


(3.20) 


+  <i(s 


1  |  eOTP 


+  spT0)  -  u 


T  (pa) 


)  |T  +  4Q 


X (pa)y 


y  x  ■ 


The  symmetry  of  Gpa  and  the  second  and  third  terms  of  the 
right  hand  side  in  p,a  gives  at  once  the  six  equat ions  of 
motion  for  spin : 


er(a)  (  6L/6  e (a)  )  =  0  .  (3.21  ) 

L  P  a  J 

The  gravitational  field  equations  are  the  remaining  ten 
equations  of  (3.20).  As  a  first  step  towards  reducing  these 
to  a  familiar  form,  the  identity  (2.16)  is  used  to  replace 

( 6 L / 6 e ^ ) e ^ a ^ p  with  terms  having  a  direct  physical  inter- 

a 

pretation.  Let  R  and  0A  denote  Ra0  .  and  0  Bt<S.  In  the 

A  3y<5  a 

present  context  of  (3.10)  and  (3.11),  the  identity  (2.16) 
reads  1 

(  sL/s  e  (a)  )  e  (a)  =  UTa  |  +  t  CT  +  QA  ( i  B  )  aRi3  (3.22) 

a  p  p  |  t  p  ApB 

modulo  the  nongrav i tat i oral  field  equations : 


^he  last  term  of  (3.22)  has  the  explicit  form 

QA ( i A B  )  ctRr  =  Qn^6RCTR  x  -  30 

ApB  P  By  o  «  py  5 
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6L/fi<(,  a  =  0  .  (3.23) 

A 

Substitution  of  (3.22)  into  (3.20)  gives  the  following  form 
of  the  grav it  at ional  field  equations : 


(  8tt  )  -  ^  /-"g  Gp  0  ~  /-~g  Tp  0  = 


tpa  +  I(s0Tp  +SpT0  +S0pT),  +  40 


(3.24) 

X  (pcr)y  +  0A(I  B)  paR 


y  X 


'A 


B 


3.4  Matter  and  Field  Decomposition  of  Energy -Momentum  and 
Spin 

Equation  (3.24)  identifies  Tpa  as  the  "correct"  (sym¬ 
metric,  covariantly  constant)  total  energy -moment urn  tensor. 
It  is  expressed  in  terms  of  a  canonical  energy  tensor  den¬ 
sity  t  a,  a  spin  flux  SpCTT  and  gravitational  quadrupole 

p 

terms.  The  presence  of  interaction  terms  in  L  implies  that 

t  0  and  Spax  will  not  in  general  be  simply  a  sum  of  free 

p 

material  and  field  parts,  but  will  contain  (respectively) 
interaction  momenta  and  spin.  As  Israel  [6a]  points  out  for 
the  electromagnetic  case: 

"This  expression  for  the  total  energy  tensor  is  the 
fundamental  result,  and  questions  about  which  part 
should  be  called  the  'electromagnetic  energy  tensor' 
are  mere  semantics  and  to  a  large  extent,  super¬ 
fluous.  However,  if  a  prescription  is  desired,  even 
though  it  be  arbitrary,  the  least  that  one  should 


. 
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demand  is  that  it  be  simple,  natural,  general,  and 
unambiguous . " 

In  the  context  of  Lagrangian  dynamics  one  immediately  has  a 
simple,  natural,  general  and  unambiguous  split  into  matter 
and  field  parts  from  the  following  considerations. 

Consider 


t 


=  L5 


-  (3  am) 

p 


aL 


3(3  am) 
a 


-  e 


(a) 


3L 


3L 


a 


3  ( e  ( f } ) 

a  cr 


A  p  3  <J> 


A  a 


If  a  convention  must  be  adopted  for  a  split  then  the  obvious 
one  is  to  designate  the  second  and  third  terms  as  belonging 
to  the  material  since  their  form  involves  differentiating 
and  multiplying  by  material  derivatives  3aam  and  e^g. 
Similarly  the  last  term,  obtained  by  differentiation  and 
multiplication  by  <f>  .  ,  is  regarded  as  field  energy-momen- 

Ci 

turn.1  The  splitting  of  t^a  is  completed  by  decomposing  L 
into  a  sum 


L  =  L1(V  't'A|a)  +  U2(V  ^  A  |  a  1  (3'25) 

in  which  =  /rgL1  and  L0  =  /^L2  represent  matter  and 
field  parts  respectively.  We  thus  have 


1  This  is  a  very  general  procedure.  If  a  set  of  variables  ^ 
on  which  a  Lagrangian  density  L  depends  are  of  two  types,  ^ 
=  (  0  A ,  nA)  ,  then  our  convention  immediately  splits  ^  A  |  p  L 
into  energy -moment urn  eA]p8L/90A]a  of  one  and  energy- 

-momentum  pa| p 3L/ 3n A| a  of  type  two. 
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=  t 


where 


p(mat) 


■,  <$ a 
1  P 


a  +  t  a 

P  (  4>  )  p  (mat) 

(3.26) 

(aam)  aL  -  aL 

p  „  a  p  /  _  \ 

(3.27) 

3(3  a  ) 
a 


3  ( e v ,  y ) 
a  a 


defines  the  material  energy -momentum  tensor  density  and 


t  ^  =  /-g  t  v  =  L  6°  -  (p  |  3  L  /  3  .  i 

P  ( )  P  ( )  2  P  A I  p  A  |  c 


(3.28) 


defines  the  canonical  energy  tensor  density  for  the  fields 

<t>  .  The  particular  split  of  L  into  L.  and  L„  is  of  minor 

A  1  1 

importance.  Equations  (3.27)  and  (3.28)  involve  the  total  L 
except  for  the  "diagonal"  terms  1^6°  and  Changing  the 

split  of  L  wi 1 1  merely  redistribute  energy-momentum  between 
these  diagonal  terms. 

The  conditions  uaa  am  =  0,  u  ua  =  -1  determine  the 


a 


a 


a 


four -ve loci ty  u^  as  a  function  of  a  am  (m  =  1,2,3)  (cf.  Sec. 


a 


2).  If  u  appears  in  L  it  is  to  be  regarded  as  such.  Noting 
that  the  second  term  on  the  right  hand  side  of  (3.27)  is 

C i 

orthogonal  to  u  only  in  its  first  index,  we  write  it  as  a 
sum  of  a  convective  four -momentum  flux  and  a  stress  term  by 
projection  (on  the  second  index)  parallel  and  orthogonal  to 
ua.  This  can  be  expressed  neatly  by  introducing  a  new 
Lagrangian  that  includes  ua  among  its  variables.  In  terms  of 
the  projection  operator 


A3  =  S3  +  U  UB 


a 


a 


a 


(3.29) 


- 
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and  L  ( am,  3aam,  e<a)  ,  e<j°B ,  *A>  *A|B,  R“By5)  we  define 


L' ( u“  ,  am ,  3  am , . 

a 


\  _  i  r  _  in  .  $  „  \ 

/  -  L\ a  ,  Aa8ga  , . . . ) 


(3.30) 


It  then  follows  that 


=  — — —  u  (3  am)  ,  — 3L'  -  -3L  act  ,  (3*31) 

8tlP  8  ( 8^  am)  a  p  8 ( 3gam)  8 ( 8a  am)  a 

giving 

-  (8  am) — — —  =  3L  uq  -  (3  am)__?_kl_  .  (3.32) 

p  8(8  am)  8  up  p  8(8  am) 

a  a 

From  (3.31)  L'  satisfies 


8L'  a 
- u 

8 


=  0 


u 

a 


8L' 


8(8  am) 
a 


=  0 


(3.33) 


It  seems  most  simple  to  introduce  ua  in  (3.30),  although  we 
could  have  considered  a  Lagrangian  dependent  on  both  3aam 
and  ua  at  the  outset  in  (3.10)  and  (3.11).  Such  a  Lagrangian 
is  largely  arbitrary  as  a  function  of  3aam  and  ua  due  to  the 

constraints  u  ua  =  -1,  ua8  am  =  0,  so  that  a  class  of 

a  a 

Lagrangians  is  associated  with  any  given  Lagrangian.  One  can 
show  that  all  members  of  a  given  class  of  Lagrangians  give 
the  same  energy -momentum  on  variation  of  the  metric  (Appen¬ 
dix  2).  The  neatest  presentation  of  the  results  is  in  terms 
of  L'  whose  dependence  on  ua  and  8  am  has  been  delimited  to 

a 

a  unique  form  by  conditions  (3.33). 

The  dependence  of  the  Lagrangian  on  the  angular 
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velocities  waD  =  e^a;a  &°{  ue  determines  the  internal  spin. 


ab  _  e(a)a  e(b) 

a|  $ 

(This  is  demonstrated  in  Appendix  10  for  a  single  particle, 


and  hence  for  the  material  as  a  whole.)  We  therefore  assume 


that  e 


(a) 

Is 


a 


equivalently,  via  e 


appears  in  the  Lagrangian  via  wab  or, 

•  (a)  -  ^(a)  ,,3 


a 


=  e 


a 


|  3  U  • 


Wi th  L' (u 


a 


(a) 


’  6  a  |V - )  =  L  ’•  *  ’ 


(a) 


a 


(3.34) 


it  fol lows  that 


-  e 


(a)  3L 


a  I  P  (a) 
3e  |  ' 
a  a 


Ci  e(i°  3L  Ii°  = 

a  I P  •  (a) 


3  e 


a 


— 

3Uf 


3L 


/ 1 


3  UPJ 


.a 


(3.35) 


(3.32)  and  (3.35)  give 


-  (3  -  e  ,  ,  , 

P  a!P  ^  tSV 


3  (3  a  ) 
a 


(a)  =  /Z 

(  ' 

a  e 


3  L  ,,o  3  L 


a  a 


1 


3  U 


-  uu  -  (3  a“> 

p  p  3(3  am) 

a 


Inserting  this  into  (3.27)  yields 


’p(mat) 


3  L 

3  U  1 


-  -  L  U 

P  1  P 


u°  + 


a  .  _  m  3  L 
LA  -  (3  a  ) - 

i  P  p 


m. 


(3.36) 


3(3aa  )j 


where  t 


p  (mat) 


=  /=it 


p (mat) 


.  Accordingly  we  define  the 


canonical  material  four-momentum  P  and  the  pressure  tensor 


P  °  as 


P  =  3  L/  3  up  -  L  u 

p  1  p 


P  a  =  L  act 
p  1  p 


O  am)..  *L.— 
p  30  am) 


(3.37) 
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so  that 


t 


a 

p(mat) 


P  ua  +  P  a 

P  P 


(3.38) 


From  the  same  considerations  as  the  first  footnote  of 
this  section  the  total  spin  flux  can  be  decomposed  into 
matter  and  field  parts.  (2.14)  and  (3.18)  give 

SP0T  =  2UX[PO]  =  2lf(I  ?)[ap]* 

A  B 


=  ?  dL  (-^[pga]g)p^a)  +2^= — (IA~)  L^Jcj> 


a  e 


(a)  a 


3cf> 


A  x 


B  [op  ] 

A  )  ^  B 


a  x 


In  terms  of  the  material  spin  S 


pa 


spo  E  2  (a) [p  3L 

(a) 

a  ] 


3  e 


and  field  spin  flux  S 


pax 

(4>) 


pax  _  2  3  L  ^  B )[  a  p  ] 


(1)) 


A 


B 


we  have 


(-g )'"2  SpaT  =  sp0  UT  +  SPPX 


(3.39) 


(3.40) 

(3.41 ) 


This  decomposition  is  independent  of  the  split  of  L  since 
(3.39)  and  (3.40)  are  in  terms  of  the  total  L. 


' 


36 


3.5  Balance  Laws 

The  Einstein  tensor  satisfies  GagiD  =  0  and  G^a^  =  0. 

I  p 

The  Einstein  equations  (3.24)  are  therefore  inconsistent 
unless  Ta^  satisfies  the  same  identities.  These  identities 

ct  3 

for  T  may  be  derived  from  the  action  principle  by  demand¬ 
ing  that  I  takes  an  extremal  value  for  the  actual  trans¬ 
lational  motion  and  spin  evolution  and  then  simplifying  with 
the  aid  of  (3.23).  Setting  61=0  for  "variation  of  world¬ 
lines",  variation  of  six  spin  co-ordinates  and  variation  of 
the  non-gravi tat ional  fields  will  thus  ensure  the  consist¬ 
ency  of  (3.24).  The  resulting  ten  equations  for  the  four- 
-momentum  and  spin  are  usually  referred  to  as  balance  laws, 
equations  of  motion  or  as  "local"  conservation  laws  for 
total  spin  and  total  four -momentum. 

The  spin  equations  have  been  given  already  in  (3.21). 
According  to  (3.20)  they  immediately  imply  T^ae^  =  0.  To 
express  them  as  a  balance  law,  making  use  of  (3.30),  (3.34) 
and  (3.39)  gives  the  spin  equations  (3.21)  in  the  following 
form: 


/  c  o  a  \  (a)  3  L  *(a)  3  L 

P  |o  [p  [p  g  4(a) 

J  a  ] 


(3.42) 


Identity  (2.11)  applied  to 


L(ua,  a 


m 


3  a 

u 
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.(a) 
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e(a).  O 

a  A 


(3.43) 


a 
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By  6 


(a)  3  L 


3  e 
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<a) 
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3  L 


3  e 


(a) 


)  )  gives 
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aL 

a  up 


(3  am) 

P 


aL 


3(3  am) 

a 


it  (i  5) 

34>  .  A  '  p 
A 


a 


V 


6  a  L 


(3.44) 


Comparison  of  (3.42)  and  (3.44)  gives  the  balance  law  for 
spin 


f  [  pa  ] 
(ma  t ) 


+  (aL/a$A)  dA~) 


B 


(3.45) 


The  equations  governing  the  translational  motion  are 
obtained  from  extremizat ion  of  the  action  integral  on 
"varying  the  wor Id- 1 i nes" 1 .  We  consider  a  1 -parameter  family 
of  tetrad  fields  e^  ;(x,e)  and  congruences  xy  ( a  ,  t,e).  We 
extremize  the  action  integral  I(e)  of  (3.10)  subject  to 
fixed  end-points  xy(am,t.,e)  =  xy ( am  ,  t . , 0 )  (i  =  1,2)  and 

ii 

5e^aVfie  =  0  (e^  attached  to  wor  Id- 1  i  nes ,  held  fixed  by 

a  a 

parallel  propagation).  For  the  infinitesimal  variation 
xy(am,t,e)  =  xu(am,t,0)  +  e£u(am,t)  the  accompanying 
absolute  variations  are: 


and  6  e 


5$  A 
A 

=  e$  , 

A  |  p 

* 

6 ( /-gd4  x) 

= 

r  a 

e? 

i 

|  a 

-gd4  x  , 

6  am 

=  0 

6(9  am)  = 

e  (  9 

am 

)EP, 

u 

p 

1  P  » 

6UU 

=  sAV,  Ua 

5 

6  (  dx ) -=  - 

e  £ 

a  |  3 

a 

U  U 

6dx 

P  I  a 

(a) 

a 

=  0  implying 

(3.46) 

S(4<a))  =  ee 
a 

(a)  pA 

x  R 

uV  +  Se“ 

ayv  6t 

O 

6 

d  0 
XFj 

the 

appl i cat  ion  of 

this 

technique 

to 

non-spinning 

,  cf 

.  [28]. 
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The  action  principle  then  gives  for  arbitrary  £y(am,  t) 


o  =  <LL  =  — 
u  de  de 


L (  u“  , am ,  3  am,e(a)  ,e(a)  ,«  )  d4x 
y  a  a  A 


Ua,  +  It  aV,  ua 

I  a  3UU  P  I  a 


- -  (3pa  i 

a(3  a  )  I 


(3.47) 


a 


+  ~v;(ej(a)  RX  uV+  'e(ah  i  11%")  +  It  t  ,  CP 
ag(a)  A  ayv  a  p  |  a  '  3$  A  |  p 

a  „ 


/^F  d  x  . 


The  first  of  (3.33),  together  with  (3.34),  implies 


Hence 


At  u“  + 

3Ua 


9  L 

3e(a) 

a 


0 


It  4K5P 

3UP  » 


a  ,  3L  *(a)  _ 
u  +  — T-  \  e  5 
a  „  ‘(a)  a  p 


3e 


a 


P  a 

u  u 


a 


(3.48) 

(3.49) 


Substitution  of  (3.49)  into  (3.47)  and  integration  by  parts 
gives  the  following  four-momentum  balance  law 


t 


a 

P  (mat)| 


a 


_  1r  s^Ya  —  $  | 

2  pagy  (mat)  +  3$  ap 

A 


(3.50) 


3.6  Einstein-Lorentz  Theory  for  Dielectrics 

To  illustrate  the  results  of  this  chapter  we  now 
consider  a  charged  dipolar  medium  and  its  interaction  with  a 
Maxwel 1 -E i nstei n  field,  described  by  a  vector  potential  Aa 
and  the  metric  tensor. 
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The  phenomenological  current1  Ja  and  the  (skew- 
-symmetric)  displacement  tensor  Ha^  are  defined  by 

Ja  e  3L/3  A  ,  Ha3  e  -4tt8L/3AoI  =  -  8  tt  3  L  /  3  F  ,  (3.51) 

a  3  I  a  ap 

in  which  it  has  been  assumed  that  A  .  appears  in  the 

a  I  3 

Lagrangian  only  through  the  skew-symmetr i c  combination  F  = 
2 3  [  ctA  6  ]  *  The  electr°magnet  /c  field  equations  are  therefore 
( cf .  eq . ( 3 . 23 ) ) 

3.  F„  =  0  ,  Ha6,  =  4ttJ“  ,  (3.52) 

[a  By]  I  3 

and  imply  conservation  of  free  charge 


Ja,  =  0  .  (3.53) 

I  a 

In  order  that  equations  (3.52)  reduce  to  those  of 
Maxwell  in  the  absence  of  matter,  L  must  reduce  to  the  free- 
-field  electromagnetic  Lagrangian2 

L  =  -(16tt)_1F  Fyv  .  (3.54) 

0  y  v 

Defining  the  electromagnetic  polarization  tensor 


MoB  e  2 3 ( L  -  L  )/3F  . 

0  a  3 


(3.55) 


^he  microscopic  current  is  ja  =  6(L-L0)/ 5A 
2  Lovelock  [27  b,c]  has  shown  that  the  most 
Ln(A  3  3  A  a ,  ga  3  )  fo  r  which  6L0/6Aa  =  Fa3>3 
(3.54)  plus  a  trivial  divergence  term.  1 


a 


=  Ja+Ma3i3 

general  ' 
is  given  by 
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leads  at  once  to  the  usual  Lorentz  polarization  relations 

HaiS  =  paB  -  47rMaB  .  (3.56) 

Since  only  in  non-conducting  material  can  dissipative 
processes  be  expected  to  be  absent,  we  therefore  demand  that 
L  satisfies 


Ja  =  3  L / 3  A  =  eN 

a 


a 


(3.57 


From  (3.4)  and  (3.53)  it  follows  that  the  charge  per  par¬ 
ticle  e  satisfies  de/dx  =  0.  Assumption  (3.57)  amounts  to 
assuming  that  the  (undifferentiated)  potential  A  appears  in 
L  only  through  a  bilinear  interaction  term 


a 


-  ^-g  eA^N  -  eN(am)Aae  5 ( 8 ^ a1 )  ( 9^a2 ) ( 3 a3 )  .  (3.58 


From  (3.58)  it  is  apparent  that  is  independent  of  the 


metric  and  therefore  does  not  contribute  to  Tpa.  Further¬ 
more,  since  L  is  the  only  gauge-dependent  part  of  L, 

(with  A  |  entering  via  F  Q,)  Tpa  must  be  gauge  invariant. 

a  |  3  a  p 

To  ensure  that  the  decomposition  of  Tpa  is  into  gauge 
invariant  parts,  some  adjustments  must  be  made  to  the 
definitions  of  matter  and  field  energy -moment urn  and  spin. 

The  redefinitions  are  simplest  when  L  is  assumed  to  be 

\  A  ) 

part  of  L  .  Then  it  is  easily  seen  from  (3.37)  that  both  P 
K  l  p 

and  a  "Kinetic  four -momentum"  p 


' 
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p  e  P  -  en A  (3.59) 

p  p  p 


are  gauge  invariant  (Appendix  3).  From  (3.28)  and  (3.40)  we 
obtain 


a 


'p  (4> ) 


=  L  2  6 

2  P 


-  ( 4  7T )  1 A  i  H 


a  a 


a 


4ttS 


p  CTT 

(40 


2A[pHo1t 


(3.60) 


so  that 


’  p  ( 4> ) 


+ 


1 


-( S  Ta 
2^  p(40 


+  S 


a  t 


( 4> )  P 


+  S 


a  t 


P  (4> ) 


T  a,  ,  -  A  Ja  (3.61) 

P ( em)  p 


where  we  have  defined  a  gauge- i nvar i ant  electromagnetic 
energy  tensor  by 


T  a 

p  (  em) 


(  4  7T  ) 


-1 


F  H 

p  a 


a  a 


+  L  5a 

2  P 


(3.62) 


From  (3.26),  (3.38),  (3.59)  and  (3.61)  we  obtain 


(-g) 


2-£  a 

P 


+  — ( S  Ta 

p  (1)) 


+  S 


a  x 


(40  P 


+  S 


a 


p  (<)0 


(3.63) 


O  |  O 

=  D  U  +  P 

'  p  p 


+  T 


a 


P  (  em) 


This,  together  with  (3.41),  gives  the  grav i tat ional  field 
equations  (3.24)  in  the  form 


( 8 tt )  1Gpo  =  Tpa  E  Tpa  +  Tpa 

( ma t )  ( em) 


(3.64) 


1,-PTO  .  ,olp  .  „opT  ,  ,_X(po)ll  ,  .A,,.  B.  pa 


+  +  S(mat)  +  s(mat))  |t  +  4Q 


yX 


+  Q  (i*  >  R 


B 
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in  which  the  tensors 


Tpa 

(mat) 


>pua  +  P 


pa 


rp  ax 
(mat ) 


SpauT 


(3.65) 


represent  the  material  fluxes  of  four -momentum  and  spin 
angular  momentum.  Their  divergences,  according  to  (3.45)  and 
(3.50),  may  be  expressed  as 


lrpax 

2  (mat )  |  T 

T  [P°  1 
(mat) 

_  F[p  Ma]a  - 

a 

4  d [ P  na]a$Y 

aBy^ 

(3.66) 

(T,  N  +  T,  , 

(ma  t )  ( em) 

)  a,  = 

P  1  o 

_Id  c$ ya 

2  pa  By  (mat) 

+  Q  , 

a  Byo | p 

(3.67) 

5 

where  we  have  used  the  result 


T  0  , ,  =  -  F  Ja  -  ^MyvF  ,  (3.68) 

p ( em)  | a  pa  2  y v  |  p 

which  follows  from  (3.62),  (3.56)  and  (3.52).  As  we  remarked 
in  Section  5,  these  equations  can  be  used  to  verify  directly 
that  the  total  energy  tensor  Tpa  is  symmetric  and  conserved. 

From  (3.62)  T.pa.  depends  on  the  split  of  L  only 
through  its  second  (diagonal)  term.  A  change  in  the  decom¬ 
position  of  L  wi 1 1  merely  redistribute  terms  between  the 
diagonal  parts  of  the  the  material  and  field  energy  tensors. 
(3.62)  differs  from  proposals  of  Abraham,  and  Einstein  and 
Laub  [2]  for  the  localization  of  electromagnetic  energy  and 
momentum  in  a  dielectric  medium.  The  split  (3.64)  of  Tpa 
into  matter  and  field  may  be  called  a  generalized  Minkowski 


. 


. 
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splitting  since  (3.62)  becomes  the  tensor  proposed  by 


Minkowski  [2]  when  L.  is  chosen  to  equal  -  (  1  6tt  )  1 F  Hap  .  The 

/  a  6 

choice  (3.54)  for  U  seems  to  be  of  special  importance.  This 


jot  3 


gives  a  T 


pa 

(em) 


that  has  been  recognized  for  some  time  in 


special  relativistic  dielectric  theory  as  being  of  special 
significance  [3a,  3b,  4a],  while  Israel  and  Stewart  [6b] 
have  recently  given  persuasive  reasons  for  its  use. 

Field  equations  (3.64)  and  (3.52)  (with  (3.56))  both 
contain  divergence  terms,  namely  the  " Bel i nf ante-Rosenfe Id" 
spin  term  and  4QX  Cpa)'M  (  ^  -jn  (3.64)  and  the  polarization 


.a  3 


yA 


current  4rr M  |  ^  in  (3.52).  This  suggests  interpreting  the 
spin  and  quadrupole  terms  as  "gravitational  polarization" 
contributions  (polarization  of  energy -momentum) .  Chapter  5 
will  explore  the  concept  of  gravitational  polarization. 


3.7  Spinning  Fluids  and  Dust 

Recalling  that  equations  (3.7)  and  (3.5)  determine  the 
number  density  n(x)  as  a  function  of  am  and  a  ara,  assume 

J  a 

that  L  depends  on  3aam  only  via  n: 

L(u\  a\  3  a1”, .  .  . )  =  L  ( ua ,  am,  n,...)  .  (3.69) 

Ct  r 

From  (3  am ) ( 3n/3 ( 3aam ) )  =  nAa  (Appendix  3)  it  follows  that 

P  P 

the  pressure  tensor  (3.37)  in  terms  of  pressure  P  is 


P 

P 


P  =  L,  -  naL  /an 

1  r 


a 


PA° 

p 


9 


(3.70) 


iCi  c  ME  1  S-f  •  y>  j-j  fW 
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The  general  considerations  of  this  paper  are  therefore 
appropriate  for  the  description  of  ideal  (dissipation-free) 
spinning  fluids  whose  spin  density  is  convective,  given  by 

Sa®y  =  S°^uY.  These  are  the  " Weyssenhof f "  spinning  fluids 

(mat) 

of  [26].  (A  much  more  complicated  (and  realistic)  descrip¬ 
tion  of  spinning  fluids  may  be  obtained  in  [6d].) 

It  was  noted  in  sections  3.4  and  3.6  how  a  change  in 
the  decomposition  of  the  total  Lagrangian  will  re-distribute 
"diagonal"  energy-momentum  between  matter  and  field.  The 
definition  of  material  pressure  P  wi 1 1  therefore  depend  on 
the  particular  split  of  L,  as  may  be  seen  from  (3.70). 
Setting  P  equal  to  zero  is  therefore  only  a  meaningful 
criterion  for  "dust"  as  long  one  is  dealing  with  a  fixed 
decomposition  of  L.  Consider  the  case  where  L2  is  taken  to 

be  the  free-field  Lagrangian  L„  =  L .  ( <j>  ,  &  .  ,  e  )  (all  the 

2  0  T  a  A  |  a  a£ 

interaction  terms  in  L  allocated  to  the  material).  From 

r 

(3.70)  the  definition  P  =  0  for  dust  is  n3L1/3n  =  L1 
giving  L1  =  nLD(ua,...)  with  Ld  independent  of  n.  Lq  is  then 
the  single  particle  Lagrangian  of  (2.18)  (modulo  -u^u0, 
factors ) . 

Finally,  when  L  is  a  function  only  of  n,  L„  =  -p(n), 

F 

from  (3.37)  and  (3.70)  we  obtain  the  usual  expressions  for 

P  and  P  for  a  non-spinning  fluid  in  a  gravitational  field. 

p 

p  =  pu  ,  P=n(8p/8n)-p  .  (3.71) 

p  p 


. 


4.  HIGHER  DERIVATIVE  COUPLING 

Chapter  3  considered  Lagrangians  depending  on  variables 

and  first  covariant  derivatives  ^Al  .We  now  discuss  the 

A  |  a 

more  general  situation  where  L  is  allowed  to  depend  on 
higher  derivatives  of  the  fields.  The  notation  of  Chapter  2, 
Sec.  4  is  used  with  a(n)  denoting  a  symmetrized  set  of 
indices.  With  the  convention  that  aL/a^  has  the  same  sym¬ 
metries  as  y  ,  the  use  of  symmetrized  indices  implies  that 

the  multipole  moments  (defined  as  a  L  /  3  <*>  A  i  ,  1  have  the  same 

^  A | a (n) 

symmetries  as  in  special  relativity. 


4.1  Generalization  of  Fundamental  Identity  (2.16) 
Consider  a  scalar  density  )  where  Y A  = 

U  ,  ,  n  =  0,  1,...)  and  extend  definitions  (2.13), 

Y A | a  (n) 

(2.12),  (2.14)  as  follows: 


•  Ao^  .  .  .  an  _  |^Aa  (n)  ^ 


9  L  /  3  ip 


A  a  ( n) 


Aa(n)  _  y  (_1)m  j_Aa(n)3(m) 

’*  ~  n  ^  I  B  (m) 

m  =  0  1  ~ 


B  N  a  r  ,  Aia(n) 


uT0p  s  Op°  I  L 

w  K  n  =  0 


B  a  (  n  ) 


t  °  =  6°L  ~  I  (n+l)ij> |  ,  ,  \ \ 

P  P  n  =  0  A I  (pa  (n)  ) 


Aaa ( n) 


The  condition  (2.11)  that  L ( ^ A )  be  a  scalar  density 


(4.1  ) 

(4.2) 

(4.3) 

(4.4) 


i  s 
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. 
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3L  (1.5)  V  +  S°L  =  0 
¥?a  a  p  5  p 


(4.5) 


With  ^  =  (^|  .  n  =  0,  1,...)  this  is 

A.  A  a  (n  ) 


l  LA2(n)d 


n  =  0 


Aot(n) 


B3  (n)  v  a 

)  ip 


B  I  3 (n)  p 


+  6  L  =  0 


(4.6) 


This  generalizes  (2.15).  The  two  steps  leading  to  (2.16) 
must  now  be  generalized.  First  note  that  repeated  use  of 
(2.5)  gives 


. Ao(n)  B3  (n)  a 

L  (.  I  .  ,  s  ~  )  ip 

Aa(n)  n  r 


LA2(n)(I  B)  %  .  ,  -  nL 

A  p  B  a  (n ) 


Aaa (n- 1 ) 


P  VB|B(n) 


A  (pa(n-l ) ) 


(4.7) 

,  ( n= 1 , 2 , . .  )  . 


Next,  we  rewrite  the  first  term  on  the  right-hand  side  of 
the  above  equation  as  a  divergence  plus  a  term  not  involving 
derivatives  of  ip  .  The  simplest  method  of  achieving  this  is 
to  note  that  (4.2)  implies 


Axa(n)  _  LAa(n)  _  , Aa(n) 


T 


_ 

* 


(4.8) 


Hence 


V  iA~  (n)/,T  B)  ai 

n  =  0  A  P  B | a (n) 


r  Axa(n)  .  Aa(n),  B  s  a 

I  <U  |T  +  L*~  )(l.  )n  ^ 


n=  0 


A  ;p  ^B|S(n) 


i  Y  , Axa(n)  b  a i 

-  ( J0L*  ~  (lA  >P*B|a(n 


)  _  y  LAxg(n)  b.  a 

)  I T  n=o  *  A  p  B | a (n) x  + 
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l  L 


Aa(n)  B  a 
(**  ^ 


n  =  0 


A  'p  a (n ) 


DTa  ,  +  Ltd/)  V 

p I T  *  v  A  p  y  B 


(4.9) 


(4.6),  (4.7),  (4.9)  together  yield 


—xa  —  a  . A ,  B .  a 

U  pit  +  lp  +  l*(ia  >p  =  0 


(4.10) 


It  will  become  clear  in  the  next  section  that,  for 
variational  purposes,  equation  (4.10)  is  more  useful  when 

written  in  terms  of  UTa  and  t  a  defined  by 

P  P 


U 


TO 


P 


l  (n+l)ij;A|  (pa(n)  ) 
n  =  0  1 


LAaTa(n)  (4<11) 


a  _ 

*P  =  5PL 


^  ^A | (pa(n) ) 
n=  U  ~ 


Aaa(n) 


(4.12) 


(4.12)  is  the  natural  definition  of  a  canonical  energy 
tensor  for  Lagrangians  containing  higher  derivatives  (cf. 
[29] ,  [15,  p.122])  . 

From  (4.8),  (4.11)  we  have  YTQ^it  =  t  a  -  t  a  so  that 

P I T  P  p 


(4.13) 


Inserting  (4.13)  into  (4.10),  we  obtain  the  generalized  form 
of  i dent i ty  (2.16) 


(4.14) 
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4.2  Generalized  Action  Integral;  Variation  of  Tetrad 

As  a  generalization  of  the  action  integral  (3.10)  we 
now  consider 


I  = 


L<',A|a(n))  d'x  +  (16")'* 


Z1?  R  d4 


x 


(4.15) 


*A  =  eaa)'  ^  A  *  R%Y6) 


(4.16) 


in  which  arbitrary  symmetrized  derivatives  of  the  field 

variables  <J>A,  R01^  ^  are  now  permitted  to  appear,  but  we 

still  assume  that  second  and  higher  derivatives  of  am  and 

„(a) 


a 


are  absent. 

Under  a  variation  of  the  tetrad  field  we  have  (Appendix 


4) 


6(e)L  =  L*6(e)^A  +  uTVrax 


( di  v 


(4.17) 


with  L*  ,  UTa  defined  in  (4.2),  (4.3)  and  (4.11). 
Define  the  gravitational  multipole  moments  q 


$y 6  A (n ) 


a 


and  associated  quantities  as  special  cases  of  general 

definitions  (4.1),  (4.2): 


<ZY6i(n)  =  ^IqAi(n>  s  3L/3RA|i(n)  =  3L/3R“6Y5lA(n)  (4.18) 


3y6A(n)  _  - - nAA(n)  _  Z  m  AA(n)B(m)  1Q\ 

°a  -  -  Ay-1)  q  '  |  6  ( m )  •  (4-191 

m=  0 
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F rom  (4.16)  we  have 

L>(e)*A  =  ( 6  L/<5  e^a  ^  )  6epa  ^  +  Q^SR0^ 

giving  (4.17)  in  the  form  = 

(«L/6e^a)  )6e<a)  +  UTap6r£T  +  o/^R0^  +  (div)  .  (4.20) 

This  generalizes  (3.13),  and  a  calculation  patterned  after 
(3.15)  to  (3.20)  gives 

(8w)-1/^  Gpa  =  (6L/8e^a))e(a)p 

(4.21  ) 

+  (i.(SpTa  +  SaTp)  -  UT^pG))  +  4QX  (pa)y 

2  |  t  |  yx  • 

We  thus  recover  the  spin  equation  in  the  form  (3.21).. 


With  =  ( am ,  e^a^  ,  <f>A,  R01^  ^  )  the  fundamental  iden¬ 
tity  (4.14)  is 


(<5L/6e<a) )e(pa)  = 

(4.22) 

UTa  I  +  t  G 

+  QA(I.B  )  °R  +  LA  ,  (I  B  )  V 

A  p  B  *  (<p)  A  pyB 

P  |  T  p 

OO 

Aa(n)  _  y 

where  =  Iq 

(-l)m(3L/34>A,  ,  ,  . 

A | a (n ) p (m)  | p  (m) 

(4.23) 

Inserting  (4.22)  into  (4.21)  gives  the  grav i tat ional  field 
equat ions 


. 
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(Sn)'1/^  Gp0  =  /=l  Tpa  =  tpo  +  i(SpTa  +  S0TP  +  S 


apt 


+  40Xtpa)u  x  +  OAilAB)p0RB  ♦  LA(+)(iaV%b  . 


(4.24) 


The  $  i  w  \  represent  the  external  fields  and  their 

A  I  A  ( n ) 

derivatives.  Although  <b  may  be  chosen  so  that  their  field 


equations  take  the  form  L‘ 


C4>) 


=  0  we  still  include  the  last 


term  of  (4.24)  for  the  following  reason:  for  the  electromag¬ 
netic  field  it  is  more  convenient  to  choose  4>  A  =  (A  ,F  _) 

A  a  a  3 

(with  (b  i  s  =  F  .  for  n=1,2..)  rather  than  <}>  .  .  ,  = 

VA  X(n)  aB  X(n)  A|X(n) 

~ 

A^^^  ( n  =  0 ,  1  .  .  .  )  ;  even  though  =  0  only  for  the 

latter  choice.  When  L  depends  only  on  first  derivatives  of 
A  ,  as  in  chapter  3,  both  conventions  for  $  give  the  same 

Gt  A 

results  with  little  difference  in  ease  of  derivation. 

Chapter  4  followed  the  more  traditional  method  of  setting 
4>  A  =  Aa .  The  gauge  dependent  t  was  then  combined  with 
the  gauge  dependent  field  spin  flux  terms  in  the  usual 
manner  to  obtain  gauge  invariance.  When  L  depends  on  higher 
derivatives  of  A^  the  first  advantage  of  using  <J>A  =  (Aa,F a^) 


is  the  gauge  invariance  from  the  outset  of  t  and 

(only  step  (3.59)  is  needed  for  complete  gauge  invariance  of 
all  parts  of  Tpa).  The  second  advantage  is  that  the  multi  - 
pole  moments  3 L/  3<f>  A  |  x  ^  =  3L/ 3Fag  |  x  ^  have  the  same  sym¬ 
metries  as  the  special  relativistic  definitions  of  the  next 
chapter . 


pai 


- 
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4.3  Matter  and  Field  Decomposition 

We  continue  in  similar  fashion  to  chapter  3,  decompos¬ 
ing  L  as 

L  =  Li'^|a(n)>  +  M\|a(n)  >  •  (4'25) 


Since  WA i  ,  N  contains  only  first  derivatives  of  a 
A  |  a (n  )  J 

e(a).  i 


m 


and 


t  fol lows  that 


t°'  =  t^  +ta 

p  p  (mat )  Lp  (<f> ) 


LRaI  (p< 


where 


and 


n  =  0 


t°,  ..  =  6°L1  -  (3  am )  3  L 

p (mat )  p  1  p 


3(3aam) 


+  0  r  O’  i 

t  / i \  -  u  L_ 

p((p)  p  2 


I  ,_,.L 


nAaa (n ) 
n))U 

(4.26) 

e(f}  3  L 

a|p^Ia) 
a  |  a 

(4.27) 

A a  a (n ) 

*(9) 

(4.28) 

The  discussion  of  the  structure  of  tp(mat)  in  chapter 
3,  equations  (3.29)  to  (3.38),  applies  without  change:  in 
terms  of  L  =  (-g)  zL  and  L(u°\...,  e^a  ,...)  = 

1  1  Oi 


a",  3aa",  e<a>,  e<f>....)  =  L!am,  A^aVJ  (4.29) 

we  have 


t  p (ma  t )  =  2  t 


p (mat) 


where 


=  p  U°  +  P  0 


(4.30) 


Pp  =  IL  *  L,  up  ,  Pp°  =  pAp  -  Opam)_ii _  .  (4.31) 

»“  S(S„am) 


In  (4.26)  tp  has  split  into  matter,  field  and  curv¬ 
ature  parts.  The  quantities  YTap  ,  UTa^ ,  SpaT  also  split: 


i f  we  define 


’ 
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U 


icr 


(<t)  P 


B \  a  n  , Axa(n) 


e  ( i  )  j  l^. 

A  P  nf0  *(+)  "B|a(n)  n=Q 


♦  '  l  (n+1 )  <i  .  , -  .  .  L 


u 


,xa 


(R)P 


.pax  _ 
J(4))  = 


Aaxa(n) 

A|  (pa(n)  )  u*(4>) 

( 4 . 3a) 

-  (IaB  )  a  l  0AXS ( n ^ R  Y  (n+I ) R  ,  0AaTS(n) 

A  p  n=0  B | a(n)  n£Q  ;  A|(pa(n)) 

(4.33) 

2U(J)a]  S(R)  =  ,  Spa  e  2e(a)Cp9L/3ela) 


'O 


] 


(4.34) 


then  (4.3),  (4.11)  give 

uTOp  =  UWP  +  UU)P  -  «L/3e<a>)e<aV  (4.35) 

and 

sp0T  =  /TJ  spauT  +  sp"  -  Sp;x  .  (4.36) 

(4>)  (R) 

As  in  chapter  3  the  decomposition  of  $pQT  is  independent  of 
the  split  L  =  L:  +  L2 .  This  split  is  only  needed  in  order  to 
partition  the  term  6^L  of  tpQ  and  thus  only  affects  the 
diagonal  parts  of  (4.27)  and  (4.28). 

To  summarize,  in  equations  (4.24)  generalizing  (3.24) 
the  total  energy  tensor  Tpa  is  found  to  consist  of: 

1.  A  material  energy-momentum  tensor  tpa  .  =  Ppua  +  Ppa 

written  as  a  convective  four -momentum  flux  and  a  pres¬ 
sure  tensor . 

2.  A  canonical  energy  tensor  for  the  fields  <f>  ,  given  by 

(4.28)  as  a  certain  combination  of  field  derivatives 

^A |  X ( n )  anc^  Qualities  formed  from  the  multipole 

moments  9  L/d(pA  j  ^  ^  according  to  (4.23). 

3.  A  divergence  formed  from  the  total  spin  flux  which 


' 
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consists  of  matter,  field  and  curvature  parts. 

4.  Contributions  from  gravitational  multi  pole  moments. 

A 

5.  A  term  proportional  to  which  is  usually  zero  by 

virtue  of  the  field  equations  satisfied  by  <j>A  and  if 
non-zero  may  be  included  in  t  . 

4.4  E instein-Lorentz  Theory 

Set  <j>  =  (A  ,  F  )  in  the  previous  section  (see 

A  CL  CL  p 

discussion  that  follows  (4.24)). 

The  electromagnetic  field  equations  are  obtained  by 

variation  of  A  Keeping  g  fixed.  The  variation  in  L  is 

a  K  aB 


5L  "  lk_6Aa  +  l  9L  .5  (  F^p  |  \(n)  ) 


3Aa  "  °8Fa3|X(n) 


A a  +  l  (-1) 
3A  n=0 


n 


a 


3  L 


,9FaB|A(n)  . 


6FaB  +  ^div) 


I  A  (n) 


Noting  that  6Fag  =  2v[^6A[3]  and  that  9  L/3  Faj3 1  x  (n)  is 
antisymmetric  in  a,B,  we  replace  6Fa|3  with  -2Vg6Aa  to  obtain 


6  L  =  3J__6Aa  +  2  l  (-Dnr 
3  A  n  =  0 


a 


3  L 


3  F 


aB | A ( n ) 


6  Aa  +  (  di  v )  .  (4.37) 


| A(n)B 


Denot i ng 


,a 


=  3  L 


,aB 


3  A 


=  -  8  TT  l  (  -  1  ) 

n  =  0 


n 


a 


3  L 


[3Fo0|X(n) J 


,  (4.38) 


I  A  (n) 


from  (4.37),  (4.38)  the  variational  principle  61  =  0  gives 
the  electromagnetic  field  equations : 


. 
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HaS|g  =  4TtJa  .  (4.39) 

In  terms  of  L0  =  -  (  1 6ir )  ’  1  Fag  FaB  (cf.  (3.54))  and  L  define 
multi  pole  moment  densities 


cxB A ( n )  . 

m  '  (*>  ;  23<L  -  Lo)/3FaB|A(n)  » 


(4.40) 


pol ar izat  ion  tensor 


„a  6 


MaB  i  l  (  -  1 )n  maB~ (n)  . 

n=  0  li<») 


(4.41 ) 


and  associated  tensors 


uo 

^aB^O)  E  £  (_  1  )  n  m°tB  Y  ( in )  £(n) 
n  =  0 


X(n) 


(4.42) 


which  imply  polarization  relations 


a$  aB  aB 

H  =  F  -  4ttM 


(4.43) 


Assume  that  the  material  is  non-conducting  by  setting 
da  =  eNa  (cf.  (3.57)).  Again  this  implies  that  both  Pp°  and 
the  "Kinetic"  momentum  pp 


P 


P 


enA. 


(4.44) 


are  gauge  invariant. 

From  (4.2),  (4.23)  and  (4.43)  we  obtain 


' 
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/  \  A  /  T  B  v  o  ,  .  »  lijCxcJr-  -i  cr « 

(-g)  L*(*)(IA  >p  =  (4,r)  H  Fap  -  J  AP  ’ 


a 


P  (<J> ) 


=  «pL2  -  \rr%  I  F  J  H 

n  =  U 


a3aa (n) 


so  that 


%(*)  +  L*(»)(Ia\\  =  ^  Ue.)  '  V°  «•«> 


where  we  have  defined  a  gauge- i nvar i ant  electromagnetic 
energy -momentum  tensor  by 


T  a  =  S°L.  +  y—  F  Haa 
p  ( em)  p  ^  4tt  a  p 


00 

— Yf  jy|Ci3(JX(n) 

2  £  aB I  (p  X (n) ) 

n-U 


(4.46) 


Inserting  (4.26),  (4.30),  (4.44),  (4.45)  into  (4.24) 
gives  the  grav it  at ional  field  equations  as 


( 8  it  )  lGpo  =  Tp  0  =  Tpa  -  +  Tpa  .  +  I(SpTCT  +  SaTp  +  SCTpT), 

(mat )  (em)  2  '  |  t 

CO 

+  4QX(pa)p\  ,  +  0A ( I n B ) pa R  -  V  R  ,  p  ,  nAp2(n)  (4.47) 

y  X  A  B  ^-.AM'aCn)) 

n  =  0  1  ~ 

where  jpa  =  ppua  +  Ppa  (4.48) 

(mat)  r 

represents  the  flux  of  material  four -momentum,  Tpa  is 

(  em) 

given  by  (4.46)  and  the  total  spin  flux  SpaT  =  (-g)_^SpaT  is 
given  by  (4.36).  Tp^^,  depending  on  all  derivatives  of  F 
and  all  multipole  moments,  generalizes  the  (generalized) 
Minkowski  tensor  (3.62). 

In  terms  of  the  multipole  moments  (4.40),  (4.41), 

(4.42)  the  field  spin  flux  is  given  by 


II  8 
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SpaT  =  (_p)-^spCTT  = 
(em)  ^  g;  b(4») 


[P 


[P 


„  0 


X  a  (n) 


Ma]XTa(n)  +  r  (  F  * 

n=0  Xyl(’ 


a  (n)  ) 


M 


(4.49) 

Xyala (n) 


5.  MULTIPOLE  EXPANSION  OF  ELECTRIC  CURRENT  AND 

ENERGY -MOMENTUM 


5.1  Introduction 

We  have  seen  in  the  earlier  chapters  how  the  varia¬ 
tional  approach  may  be  used  to  derive  the  fully  covariant 
dynamics  of  spinning  polarized  media.  By  not  specifying  any 
particular  Lagrangian  we  have  created  a  framework  into  which 
any  detailed  model  must  fit.  Of  course,  the  price  we  pay  for 
this  is  the  lack  of  physical  insight  to  be  gained  from  the 

formal  definitions  such  as  P  =  3L/3va  and  Mae  = 

a 

23(L-L0)/3Fa3.  The  objectives  of  this  chapter  are  to  gain  a 
clear  understanding  of  the  meaning  of  the  various  dynamical 
quantities  such  as  pa ,  Sa&  and  MaB  appearing  in  earlier 
chapters  and  to  explore  the  concept  of  "gravitational  polar¬ 
ization"  (polarization  of  energy -moment urn)  in  detail. 

We  shall  derive  the  same  dynamical  laws,  albeit  in 
special  relativity,  by  showing  first  that  the  structure  of 
charge  and  current  for  a  composite  particle  (extended  body) 
may  be  summarized  in  certain  uniquely  defined  quantities 
called  multipole  moments.  These  are  tensors  generalizing  the 
electric  and  magnetic  dipole,  quadrupole  moments  etc.  of 
non-relativistic  electrodynamics  and  they  combine  in  a 
certain  way  to  form  the  polarization  tensor. 

Starting  from  the  microscopic  form  for  the  current  of  a 
system  of  point  charges,  we  shall  show  that  a  multipole 
expansion  results  in  the  decomposition  of  current  into  an 
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overall  flow  and  a  polarization  part.  This  is  one  of  the  two 
steps  in  the  construction  of  the  macroscopic  Maxwell  equa¬ 
tions  for  polarized  media  from  point  charge  theory.  The 
other  step  is  to  suitably  define  an  averaging  scheme  to 
obtain  smoothly  varying  quantities.  (This  thesis  explores 
only  the  multipole  formalism  and  not  the  averaging  pro¬ 
cedures  .  ) 

A  non-relativistic  derivation  of  the  macroscopic 
Maxwell  equations  using  spatial  averaging  may  be  found  in 
[31]  (reviewed  in  Jackson  [32],  see  also  references  cited 
there).  A  detailed  analysis  of  the  relativistic  formulation 
is  contained  in  de  Groot  and  Suttorp's  book  [5]  (chapter  V). 
They  achieve  smooth  averaged  quantities  by  assuming  the 
composite  particle  density  is  large  enough  that  a  smooth 
distribution  function  exists.  This  reference  also  derives 
composite  particle  equations  of  motion,  balance  laws  and 
discusses  in  detail  the  3+1  representation  of  the  various 
dynamical  quantities. 

We  extend  the  multipole  expansion  method  by  showing 
that  it  may  be  applied  to  fluxes  more  general  than  just  the 
electric  4-vector  current.  Applied  specifically  to  the 
microscopic  energy -moment urn  tensor  (4-momentum  flux),  the 
method  leads  to  a  gravitational  polarization  tensor  Na^y 
completely  analogous  to  the  electromagnetic  polarization 
tensor  Ma^.  (The  relation  between  NaBy  and  the  spin  terms  of 
the  gravitational  field  equations  will  be  discussed  in 
section  6 . ) 


■ 
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In  formulating  the  single  particle  dynamics,  defini¬ 
tions  of  four -momentum  and  spin  arise  that  make  the  equa¬ 
tions  of  motion  and  balance  equations  appear  in  their  most 
simple  form  (the  same  equations  as  those  of  the  Lagrangian 
formulation).  Finally,  the  expansion  of  the  microscopic 
energy -momentum  is  shown  to  give  the  same  decomposition  of 
total  energy -momentum  into  material,  electromagnetic  and 
polarization  parts  as  the  Lagrangian  result,  eq.(4.47). 

We  will  therefore  see,  for  media  whose  spin  is  ultim¬ 
ately  orbital  in  nature,  that  the  microscopic  (symmetric) 
total  energy -momentum  leads  to  a  macroscopic  total  (sym¬ 
metric)  energy -moment urn  tensor  of  which  the  asymmetric  four- 
-momentum  flux  is  only  a  part.  This  part  describes  the 
overall  ("gross")  flow  of  four -momentum  just  as  the  phenom¬ 
enological  electric  current  describes  the  overall  flow  of 
charge,  while  the  total  energy -momentum  and  the  total 
electric  current  also  contain  polarization  parts. 


5.2  Classical  Microscopic  Model  of  Matter 

The  microscopic  picture  is  a  cloud  of  structureless , 
charged,  point  particles  i,  each  with  rest  mass  n^,  charge 
e  ,  world- line  za(s.),  normalized  four -ve  loci  ty  ua(s.)  = 

i  ii  ii 

dz'Vds  .  and  four -momentum  p^  =  m.ua.  The  (special  relativ- 

ii  iii 

istic)  microscopic  electromagnetic  field  equations  are 


V 


a  3  _ 


-  j“(x)  =  I 


e.ua  6 4 ( x-z  .  ( s  . ) )  ds. 

ii  ii  i 


(5.1  ) 
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where  6 4 ( x )  denotes  the  four-dimensional  delta  function.  The 
symmetric  microscopic  material  energy -moment urn  tensor  (four- 
-momentum  flux)  is 


t 


a  (3 

(mat) 


(X) 


m . ua  ue  6  4 ( x-z . ( s  . ) ) 

iii  ii 


(5.2) 


If  the  point  particles  are  bunched  together  into 
distinct  stable  groups  ,  so  that  each  one  is  part  of  a 
composite  particle  K  (extended  body),  then  our  physical 
picture  is  a  medium  consisting  of  particles  k  with  spin  and 
other  structure.  Let  Z^(t  )  denote  some  choice  of  central 
(reference)  world-line  for  particle  k,  with  four-velocity  Uk 
=  dZk/dxk.  If  we  neglect  the  charge  structure  of  each 
particle  k  by  assuming  that  the  total  charge  e,  =  I  e  ■ 

R  i  of  k  1 

lies  entirely  on  Zk,  we  have  the  microscopic  current  ja(x) 
approximated  by 


J“(X)  =  I 
k 


SkUk  {4(X-Zk(Tk))  dTk 


(5.3) 


This  part  of  ja(x)  is  due  to  the  motion  as  a  whole  of  the 
composite  particles.  ja(x)  consists  of  Ja(x)  and  a  "polariz¬ 
ation"  part  that  describes  the  contribution  from  the  detail¬ 
ed  charge  structure  of  each  particle.  We  may  also  approx¬ 
imate  t^e  .(x)  with 

(mat ) 


Tct3 
( ma  t ) 


(x)  -  l 

k 


p“u®  64(x-Zk(Tk))  dtk 


(5.4) 
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in  terms  of  the  total  four -momentum  pk 


I  m.ua  of  each 


i  i 


i  of  k 


K.  The  meaning  of  the  above  expression  is  clear:  it  is  the 
flux  of  four -momentum  ( four -moment urn  current)  of  the 
composite  particles  idealized  as  point  particles.1 

5.3  Notation 

From  (5.1)  to  (5.4)  it  is  sufficient  to  consider  the 
multipole  expansion  of  the  current  and  energy -moment urn  due 
to  a  single  particle.  Many  particle  results  are  then  infer¬ 
red  simply  by  summation,  first  over  i  for  fixed  K  to  obtain 
composite  particle  quantities,  then  over  K.  We  therefore 
consider  a  single  point  particle  with  charge  e,  rest  mass  m, 
world- line  za(s),  normalized  four-velocity  u°  =  dza/ds  and 
four -momentum  mua.  We  investigate  the  particle's  "eccentric 
dynamics"  with  respect  to  some  (arbitrarily  chosen)  time- 
-like  reference  world-line  Za(x)  with  normalized  four- 


velocity  Ua  =  dZa/dx. 


a 


ii 


7  d  ,  y 

Z  (x) 


Particle 
world-1  i  ne 


Reference' 

world- 1 ine 


M5.3)  and  (5.4)  must  be  averaged  to  obtain  smoothly  varying 
spacetime  functions.  For  example,  whenever  a  smooth  distrib¬ 
ution  function  y(x,n)  depending  on  phase  space  variables  n 
may  be  defined,  averages  of  (5.3)  and  (5.4)  are  obtained  by 
replacing  the  summation  over  k  and  integration  along  Z“  with 
integration  over  phase  space: 


,a  ,  N 

<  J  ( X) > 


u(x,fi) e(n) Ua(ft)  d n 

u  (x ,n) pa (n) U3 (n)  d a 


( See  [5,6c] ) . 
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Any  monotonic  function  s  =  s(t)  defines  a  connecting  vector 


CL  *  \  CL  *  v 

a  (x)  =  Z  (t)  -  Z  (x) 


(5.5) 


(We  do  not  restrict  s(t)  with  conditions  such  as  aaU  =  0 

a 

since  the  multipole  expansion  may  be  carried  through  for  any 
central  world -line  and  any  world- tube  "slicing". 

Let  aa  denote  daa/dx.  Let 


I  (fCO)  = 

n 


f (t) (a . 3)n  6 4 (X-Z (t) )  dx 


(5.6) 


a  _  „  v  n  _  a  i  an 


where  3  g  =3g/3x  =  g  and  (a  .  3) “  =  a  1  .  .  .a  “3  ...3 

o.  & » a  a,  a 

1  n 

I  has  the  following  three  properties: 


In(cif(x)  +  C  2  g  (  T  )  )  =  C1In(f(x))  +  C2In(g(x))  ,  (5.7) 


3  I  ( f ( x ) a  )  =  I  ,  (f (t)  ) 
an  n  +  1 


(5.8) 


3  I  (f(x)Ua)  =  I  ( df ( x ) /dx )  +  n3  I  (f(x)aa)  .  (5.9) 

an  n  a  n- 1 


Proof  of  ( i  i  i  )  : 


3  I  (f(x)Ua)  =  f(x)(a.3)n((-3  /3Za)64(X-Z(x)) )  Ua  dx 

a  n 


“n 
d  x 


(f(x)a  1  .  .  .  a  )3  ...3  64(x-Z(x))  dx 


a  l  an 


=  I  (df(x)/dx)  +  n3  I  (f(x)aa) 

a  n-1 


n 
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According  to  (5.9), 


3  I  (f(T)o°US)  =  I  (fa°) 
P  n  n 


+  I  (faa)  +  n3  I  (f aaa^) 
n  p  n  - 1 


where  the  dot  denotes  differentiation  w.r.to  t. 
Using  (5.8)  and  summing  gives 


°°  ,  ,  N  n  +  1  0  00  /  i  n  n+1 

ln Tsnr 

n  =  (J  n=  U 


-  In  (fa01)  +  8 


I 


(-D 


n  +  1 


a  3 


n  =  1 


(n+1) ! 


I  ( f (a  a 
n  —  l 


,  a  •  3  .  \ 
+  na  a  )  ) 


or 


o  =  -  i„(foa)  + 


l 

n  =  0 


(-D 


n  +  1 


(n+1) ! 


I  (fa 
n 


a 


(5.10) 


+  3 


(-1) 


n 


Jo  (n+1)! 


I  (f(aaUB 

n 


1  • a  3  ,  (n  +  1)  a  •  3S  \ 

+  (^T2)a  a  +  ( n+2*)  °  CT  )} 


5.4  Multipole  Expansion  of  Current  and  Energy -Momentum 
Having  developed  some  notation,  we  now  proceed  to 
expand  the  microscopic  fluxes  in  powers  of  oa  to  produce  a 
splitting  into  "gross"  and  polarization  parts. 

The  current  and  energy -momentum  of  the  particle  are 
given  by 


.  a 
J 


(X) 


eua(s)  64(x-z(s))  ds 


5 


(5.11) 
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muau^  64(x-z(s))  ds 


These  are  examples  of  a  general  flux 


fa(x)  =  f(s)ua(s)  s4(x-z(s))  ds 


From  za  =  Za+aa  we  have 


a,  ,  a  .  ,  a  .  ,  .  ,  ,  i  iOt  .  •  a  N  . 

u  ds  =  dz  =  (dz  /dT)dx  =  (U  +o  )dx 


and 


Inserting  (5.14)  into  (5.13)  gives 


fa(x)  = 


f(T) (Ua(x)+aa(T) )  64(x-z(t))  d 


(5.12) 


(5. 13) 


(5.14) 

(5.15) 


Expanding  6 4 ( x - z)  =  s4(x-Z-a)  as 

CO  ft 

64(X-Z(T))  =  l  -^p-  (a  .  8)ns4  (X-Z(T)  )  (5.16) 

n=  0  n' 

gives,  in  the  notation  of  (5.6) 

00  n 

fa(x)  =  l  1  (f(^)  (Ua+aa)  ) 

n  =  0  n'  n 

Let  Fa  denote  the  lowest  order  (a0)  term  in  the  above: 


Fa  ( x)  =  I  o  (  f  Ua )  = 


f (x)  6 4 (x-Z  (t)  )  dx 
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Then ,  with  (5.8), 


fa(x)  =  Fa ( x)  +  I0(faa)  +  l 


(-1) 


n  +  1 


n=  0 


(n  +  l)  !  36In(f(U“^)aB) 


( 5 . 10)  then  implies 


fa(x)  =  Fa ( X)  +  3  XaB 


(-1) 


+  >  tsfiTT 

n  =  (j 


where 


X“B  ■  l  CnTT)T  In(2fo[a(UB1+^>;B1)) 


n=  0 


(n  +  2) 


a 


Setting  f  =  e  and  f  =  mu  respectively  gives 


.  a  ,  .  ,  a ,  N  „  ..a  3 

J  (X)  =  J  (X)  +  3  ^  M 


(5.17) 


tf  ,  (x)  =  Ta6(x)  +  3  Na0Y  +  flaB 
(mat)  y 


(5.18) 


wi  th 


Ja(x)  =  I  n  (  e  U  a  )  = 


ella  64(x-Z(t))  dr 


(5.19) 


T“B(x)  =  I0(muaU6)  = 


muaUB  s4(x-Z(t))  dx  ,  (5.20) 


.a  3 


n 


t  (  orx  ,  ( n  + 1 )l  $  1  •)  ) 


e  I  (Sfr  In(2a-(U-  + 


n=  0 


(n  +  2) 


(5.21 ) 


N“BY(x)  ■  l  -^ifr  In(2muaolP(UYJ+^^o’,J)) 

n=  0 


n 


a  [  3  ,  ,,y  ]  ,  (n  +  l)  •  y  ] 


(5.22) 


and 


- 
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and 


A 


a  3 


I 


(-D 


n  +  1 


n=  0 


(n+1) !  ^  n  ^ 


d (mua)  3 


dt 


-  O  p  ) 


(5.23) 


ct  8  •  i 

M  is  the  polarization  tensor.  It  is  constructed  from  the 
particle  multi  pole  moments 


a3A (n)  ,  .  _ 

m  ~  ( t)  - 


2  e 


(n+1) ! 


[a  /  ,.3  ]  (n  +  1)  -  3]  \  X  (n) 

+  (n+2)  IO~ 


(5.24) 


and  associated  multipole  moment  densities 


ma|+(n)(x)  =  I0(maB^(n>(T))  = 


m 


aB^(n)(x)64(x-Z(T))dT  (5.25) 


according  to 


M 


a  3 


l  <-Dn3 

n-  0 


A  (n) 


maB^(n)(x) 


(5.26) 


which  is  equivalent  to  (5.21). 

(5.26)  is  the  special  relativistic  form  of  (4.41).  (5.24) 
and  (5.25)  therefore  give  a  clear  physical  meaning  to  (4.40) 
in  terms  of  the  detailed  particle  structure. 

A  in  (5.18)  arises  from  the  integration  by  parts  in 
identity  (5.9).  It  has  no  counterpart  in  (5.17)  since  de/dx 
=  0.  If  the  particles  are  in  free  motion  so  that  AaB  van¬ 


ishes,  then  (5.17)  and  (5.18)  are  the  decomposition  of  ja 
and  into  "gross"  and  "polarization"  currents.  In 

general  when  matter  and  field  interact  we  have  d(mua)/dx  and 
Aa^  non  zero.  Aal3  is  a  field-matter  interaction  term  that 


itself  must  be  split  three  ways  into  material,  field  and 


■ 
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polarization  parts. 

5.5  Spinning  Multipole  in  an  External  Electromagnetic  Field, 
Balance  Laws. 

The  Lorentz  force  equation  is 


d(mua)/dx  =  efaD (z) dze/dx 

p 


(5.27) 


In  the  expansion  (5.17)  the  multipole  moments  (5.24)  appear 
as  the  basic  entities  summarizing  the  four  current  "struc¬ 
ture".  It  is  possible  to  express  the  above  Lorentz  force  in 
terms  of  these  and  derivatives  of  fag  at  Za ,  together  with  a 
total  time  derivative  along  the  reference  world- line.  This 
can  be  seen  from  the  following  identity  (Appendix  5)  valid 
for  a  funct ion  f ( z ) : 


f(z)dz^/dx  = 


(5.28) 


Setting  f(z)  =  efa^(z)  then  gives  the  multipole  expansion  of 
(5.27)  as 


(5.29) 


where 


. 


68 


a  a 

p  Emu 


-  l 

n=  0 


-a 


(n+1) ! 


(Z) 


X  (n) 


3  X  ( n) 
a  a  ~ 


(5.30) 


ihe  total  time  derivative  of  (5.28)  has  been  combined  with 

a 

mu  in  (5.29).  This  gives  the  equations  in  their  simplest 
form,  with  the  force  determined  from  the  fundamental  quan¬ 
tities  Ua  and  of  the  previous  section.  The  resulting 

equation  is  the  special  relativistic  limit  of  the  Lagrangian 
equation  (2.29) . 

To  obtain  spin  equations  of  motion  consider 


=  i[ctmu8]  +  o[a^(>r,u81  ) 


=  mi3aU8^  +  ea^af8^  (z)(dzT/di) 

y 


(5.31 ) 


from  (5.15)  and  (5.27).  The  multipole  expansion  of  the 
second  term  in  the  above  may  be  found  in  Appendix  6.  With 
spin  angular  momentum  defined  as 


o°  ,  x 

Sa8  =  2a[a(muB1  -  l  f1 3 1  (Z)  ,,  ,aYa^(n))  (5.32) 

f:  (n  +  2)  !  y  » X  (n) 

n=  U 


the  expansion  of  (5.31)  takes  the  form 


■|dSaB/dT  =  p[aUBl 


(5.33) 


l  f[a  (Z), 


n=  0 


y 


X  (n) 


m 


3 ] y X  (  n) 


00  [  *  •  1 
+  l  (n+1 ) f a  (Z)  _  ,  ,mY5B~(n) 


n=  0 


Y  »  6  X  (  n  ) 


which  is  the  special  relativistic  limit  of  (2.30).  In 
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similar  fashion  to  the  definition  (5.30)  of  four -momentum, 

some  interaction  angular  momentum  that  appeared  in  the  spin 

equations  as  a  total  time  derivative  has  been  designated  as 

"material".  The  couple  acting  is  then  determined  from  U  and 

ma6X(n)  ^or  a  g-jven  electromagnetic  field. 

Equations  (5.29)  and  (5.33)  may  be  cast  in  the  form  of 

balance  laws.  Define  a  material  energy-momentum  tensor 
Tag 
(mat)  ' 


ja£ 

(  ma  t ) 


(X) 


MpV)  = 


.a  ,,B 


ir  64(x-z(t))  dx 


(5.34) 


(5.9)  and  (5.29)  immediately  give 


a  Ta  ^ 
d  3  (mat) 


f“  (x)  JB(x) 

P 


°°  .  . 

I  f“0(x)  w  . ~ ( n) ( x) .  (5.35) 

Zn  &  >yX (n) 


The  right  hand  side  may  be  written  as  a  divergence  (Appendix 
7 ) ,  giving 


8 


3 


(T 


a  3 

(mat) 


+ 


T 


a  3  \ 
(  em) 


0 


(5.36) 


where 


■  a  3 

(  em) 


x  =  ( 4  tt)  1(fa  fBy  -  ygaefy6f  )  -  f a  M 

y  4  y  o 


a  |U$Y 

Y 


CO 

-  ±  l  f  ,<x)  “MX«X(n)B  .  (5.37) 

2  ^  _  yo  »  A (n) 
n=  0  ' 


(5.37)  is  the  special  relativistic  limit  of  (4.46)  where  L 


•  Lo- 


■ 
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The  angular  momentum  balance  law  is  obtained  in  the 
same  way.  Let 


(mat) 


i0(s“V) 


$aBUY  s4(x-Z(t))  dx  .  (5.38) 


Then  from  (5.9)  and  (5.33) 


T  -  c ct  3 y  _  T  *  /  c  ^  ^  \  t  [  3  ] 

23YS(mat)  ~  2'o(dS  /dx)  =  T(„t) 


-  I  f[“  (X) 


n=  0 


y  ’X(n) 


m 


3 ] yX ( n) 


[ 


(x)  +  X  _  (n+l)  fCtT  Cx)  _ 


n=  0 


6  X  (  n) 


y  6  X ( n) 3  )  N 
m  '  ~  ( X) 


The  last  two  terms  are  expressible  as  the  sum  of  T;a^  and  a 


(  em) 


divergence  (Appendix  8).  In  terms  of 


<-a3y 
(  em) 


(5.39) 


l 

n=  0 


f 


[  a 


k  >  X  (  n  ) 


KB]Kj(n)Y  _  (n  +  1)| 


a 


k  j  y  X  (  n  ) 


M 


k  y  X  (  n  )  y  3 


1 


we  have 


1 

2 


8 

Y 


(S 


a  3  y 

(mat ) 


+ 


c-a3y  x 
(em)  ; 


(T 


(mat) 


T ( em)  ^ 


[  ot  3  ] 


(5.40) 


(5.39)  is  the  special  relativistic  form  of  (4.49). 
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5.6  Localization  of  Energy  Momentum. 

Having  obtained  the  same  equations  of  motion  and 
balance  laws  as  in  the  Lagrangian  theory,  we  now  consider 
the  derivation  of  energy -momentum  localization  as  embodied 
in  the  gravitational  field  equations  (4.47).  The  microscopic 
total  energy  tensor 


(  mat) 


+ 


.  a3 
(em) 


(5.41 ) 


is  given  by  summing  (5.12)  and  the  microscopic  electromag¬ 
netic  energy -momentum  tensor 


tf  .  (X)  =  (4.)-1{f^fBY  -  igoBfYSfY 4)  •  <5-42> 

All  three  tensors  in  (5.41)  are  symmetric.  We  wish  to  show 
that  the  multipole  expansion  of  (5.41)  is  the  special  rela¬ 
tivistic  form  of  (4.47). 

The  total  energy-momentum  tensor  is  often  defined  from 
the  balance  equations.  Any  conserved  tensor,  for  example 
T°^  +  T°t^  .  is  then  a  candidate.  The  balance  laws  merely 

(mat)  (em) 

state  however  that  certain  tensors  have  vanishing  diver¬ 
gence.  To  these  an  arbitrary  curl  may  be  added  without 
affecting  the  balance  law,  so  the  balance  laws  do  not 
provide  a  unique  prescription  for  localization. 

At  the  microscopic  level  equation  (5.41)  uniquely 

specifies  the  total  energy  momentum  distribution.  Only  the 

a  3 

(averaged)  multipole  expansion  of  t  may  be  legitimately 
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referred  to  as  the  true  distribution  of  total  energy-momen- 

rv  R 

turn,  which  will  be  symmetric  and  conserved  like  t  .  From 
(5.36)  and  (5.40)  we  have  that  T“6  .  +  N  is  conserved 

(mat)  (em) 

but  asymmetric,  while  one  may  easily  verify  that  the  tensor 
JaS  defined  by 


Ja3  =  (T 


,  v  +  T  /  Ja3  +  ^  (Say3  +  s6ya  +  S3ay)  ,  (5.43) 

(mat)  ( em)  2  y 


_  qOtgy  <-a3Y 

"  (mat)  (em) 


(5.44) 


is  both  conserved  and  symmetric.  This  is  the  simplest1  sym¬ 
metric  conserved  tensor  one  may  construct  from  the  balance 

laws.  Since  +  T“e  .  ,  Ja^  and  tae  are  all  conserved  they 

(mat)  (em) 

must  differ  from  each  other  by  a  curl.  We  have,  for  some 

Aa(BY)  =  0f 


.  a  3 


=  (T 


(  mat ) 


+  T(em)) 


a  3 


+ 


3  A 

Y 


a  3  y 


(5.45) 


(5.45),  (5.44)  and  (5.40)  imply  3  A[“b1y  =  -k  SaBy,  giving 

Y  2  y 


A [a3  ]  Y 


lc  a  3  Y 
2  ^ 


+  f3sQaBYS 


(5.46) 


for  some  Q“ByS  satisfying  Q(aB)y6  =  0  =  Q“B(y6). 

Any  third  rank  tensor  that  is  antisymmetric  in  its  last  two 
indices,  such  as  AaBy,  may  be  rewritten  as 


1  The  definition  of  Ja^  contains  only  those  tensors  that 
appear  in  (5.36)  and  (5.40). 


I  .  f 
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Aa6y  =  A[a3]y  +  A[ya]B  +  A[yg]a  .  (5.47) 

(5.45),  (5.47),  (5.46)  and  (5.43)  imply  that  the  required 
expansion  of  tag  takes  the  form 

ta6  =  jag  +  gy(aB)6  ( 5 . 48 ) 

>  y  5 

(5.48)  is  the  special  relativistic  limit  of  (4.47)  and 
completes  the  comparison  with  Chapter  4.  This  chapter  has 
provided  physical  interpretation  of  the  formal  definitions 
of  the  Lagrangian  theory,  albeit  in  special  relativity.  The 
fluxes  of  four -momentum  and  spin  acquire  a  direct  operation¬ 
al  significance  in  terms  of  the  distribution  of  spinning 
particles.  Localization  (4.47),  formally  defined  by  varia¬ 
tion  of  the  metric,  has  been  shown  to  be  the  (averaged) 

CL  8 

multipole  expansion  of  the  microscopic  t  .  The  present 
derivation  of  (5.48)  is  tailored  for  those  not  working  in 
general  relativity,  for  whom  the  definition  of  energy-momen¬ 
tum  as  the  variational  derivative  of  L  with  respect  to  the 
metric  may  have  little  physical  appeal.  The  concept  of 
"gravitational  polarization"  has  been  investigated  by 
comparison  with  the  traditional  account  of  electromagnetic 
polarization.  The  Bel i nf ante-Rosenfeld  spin  terms  ,  together 
with  the  gravitational  multipole  contributions,  have  been 

Ot  3 

shown  to  be  the  counterpart  of  M 


' 
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6.  MULT I POLE  ANALYSIS  IN  CURVED  SPACETIME 


6.1  Introduction 

This  chapter  will  generalize  the  multipole  analysis  of 
the  previous  chapter  to  curved  spacetime.  Trying  to  general¬ 
ize,  step  by  step,  each  equation  of  the  previous  chapter 
proves  to  be  an  extremely  arduous  task  because  the  chapter 
makes  great  use  of  the  simplicities  afforded  by  special  rel¬ 
ativity  (commuting  of  partial  derivatives,  vector  nature  of 
(Minkowski)  co-ordinate  differences,  "non-local"  character 
of  vectors).  The  following  questions  must  be  answered  before 
one  can  even  consider  a  Taylor  expansion  and  subsequent 
manipulation  of  derivatives: 

1.  how  is  the  "connecting"  vector  a  to  be  defined? 

2.  how  is  the  particle  four-momentum  (a  vector  field  along 

C ( 

the  central  world- line  Z  )  to  be  defined  in  terms  of  the 

y  y 

energy-momentum  distribution  (p  along  z  )  and  the 
applied  fields? 

In  special  relativity  the  connecting  vector  may  be 
identified  with  the  geodesic  (straight  line)  path  joining 
the  two  points.  To  generalize  this  unambiguously  to  curved 
spacetime  one  must  consider  only  points  zU ( t )  in  a  normal 
neighbourhood  of  Z  (t)  ( non- focus i ng  of  geodesics)  and  then 

CX  OL  Oi 

set  a  equal  to  al  where  1  is  the  unit  tangent  vector  at 
Z  to  the  geodesic  Zz  (of  length  o). 

In  the  definition  of  four -momentum  (and  spin)  one  needs 
a  criterion  for  selecting  the  type  of  transport  of  p^  from 


- 


75 


y  01 

z  to  Z  ,  from  all  possible  choices  that  reduce  to  parallel 
propagation  for  vanishing  curvature  [12],  As  Dixon  has 
shown1,  parallel  propagation  is  not  the  ideal  choice. 

The  approach  of  Chapter  four  does  not  generalize  easily 
to  curved  spacetime  because  their  is  no  simple  fully  covar¬ 
iant  way  of  expanding  (and  then  manipulating)  tensors.  A 
suitable  method  may  be  inferred  by  noting  that  scalars  are 
the  only  tensorial  objects  one  may  expand  simply  in  a  fully 
covariant  way,  and  the  relevant  scalar  from  which  the 
dynamics  may  be  derived  is  of  course  the  Lagrangian.  The 
first  effect  of  considering  curved  spacetime  is  therefore 
the  identification  of  "eccentric  Lagrangian  dynamics" 
(expansion  of  L)  as  the  only  conceptually  simple  approach. 
The  second  effect  is  that  qa^T^£(n)  (ancj  pa  t  s0^,  ma^T(n)) 
no  longer  needs  to  be  "synthesized"  as  in  the  previous 
chapter,  it  is  given  by  differentiation  of  the  expanded  form 
of  L.  One  merely  turns  the  Lagrangian  "crank". 

We  therefore  consider  the  eccentric,  Lagrangian 
dynamics  of  a  particle.  From  a  Lagrangian  standpoint  the 
expression  p  =  9L/9va  enters  automatically  as  the  natural 
definition  of  four -momentum.  This  avoids  comparing  the 

Oi  6 

relative  merits  of  various  proposals  for  pa  and  S  (and 
their  justifications  in  special  situations),  standard 

Lagrangian  formalism  by  itself  provides  the  answer  to 

QL  3 

question  2.  A  natural  definition  for  S  also  appears  in  the 

^ixon  defines  four -momentum  and  spin  so  that,  for  each 
syrmietry  of  spacetime  that  also  preserves  the  the  electro¬ 
magnetic  field,  a  correspondi ng  component  of  four -momentum 
or  spin  is  a  constant  of  the  motion. 


. 
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equations  to  be  derived  and  is  adopted.  These  definitions 
are  compared  with  Dixon's  in  Section  6  and  are  found  to  be 
the  same.  This  establishes  the  definitions  on  a  completely 
firm  foundation  from  two  viewpoints. 


6.2  Expansion  of  Lagrangian 

As  long  as  one  is  dealing  with  scalars,  values  at 
different  spacetime  points  may  be  compared:  hypersurface 
integrals  are  covariant  and  Taylor  expansions  may  be  written 
covar i ant  1 y 1 .  The  scalar  of  fundamental  importance  is  the 
Lagrangian.  For  a  point  particle  in  an  electromagnetic  field 
the  action  integral  is 


I 


-m 


d  s  +  e 


a.  dzx 

A 


(6.1) 


where 


dzxdzpN 
Xudt  dt 

/ 


eA  ( z) dzX 
X  dt 


(6.2) 


in  terms  of  a  scalar  parameter  t  along  world-line  z  . 

L  is  to  be  written  in  terms  of  a  reference  world-line 

o 

Za .  Consider  a  two-space  zy  ( t ,  c)  with 


VU  =  3zU(t,a)/3t  ,  ly  =  3zy(t,a)/3a 


1  For  points  z,  Z,  in  a  normal  neighbourhood  of  eachother , 
one  may  write  the  expansion  of  a  scalar  field  <f>  ( x )  as 
$ ( z )  =  y(1/n!)aStnJ+(z) |s(n)  where  aa  is  the  "geodesic" 
vector  ’joining"  Z  to  zf  a  =  ala- 
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za(t,  0)  =  Za(t)  ,  Va  =  dZa/dt  ,  (6.3) 


y  y  , .  .  ,  a  a 

u  =  v  dt/ ds  ,  U  =  V  dt/dx  , 


Ql 

(where  x  and  s  are  the  proper-times  along  Z  and 
zy ( t ,  a  =  const ) ) . 


Initially  zy ( t ,  a)  with  t  constant  are  an  arbitrary  set  of 
spacelike  curves  joining  Z  (t)  and  z  (t,  a)  with  a  any 
parameter  along  them.  Indices  a,  3,  y,...  will  be  used  for 

Ci  X 

tensors  at  Z  (t)  to  distinguish  them  from  tensors  at  z 

a 

which  will  have  indices  X,  y,  v,...  .  Z  (t)  represents  the 
reference  world- line  and  z^ ( t ,  a)  wi th  a  constant  the  actual 
particle  world-line.  From  (6.3)  we  have 


6vX/6o  =  SlX/6t 


(6.4) 
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To  expand  (6.1)  we  need  to  write  the  scalars 


¥(a)  = 


(“9yVvUvV)^  dt 


(6.5) 


$(a)  = 


Ay( z( t ,a) ) vP  dt 


(6.6) 


,a 


in  terms  of  tensors  at  Z 

Consider  first  the  expansion  of  the  scalar  iM  t ,  cr ) 

where 


^  ( t ,  a )  =  -(ds/dtf  =  VyV^(t,a)  ,  'F(ct)  =  (-ip)  2  dt  .  (6.7) 


Making  use  of  (6.4), 


dip/d  o  =  2vy6vy/6a  =  2vu6ly/6t 


y 


(6.8) 


2,„  X.,  Vl,62iy 


1  °  ^  6  v  5  1 

2  3a2  6a  6  t 


6a6  t 


61y  6 1 y  +  vlVR  y  rvXl^  +  v..  6  2 1  ^ 

M  V  A  Cj  M 


6  t  6  t 


6  t6a 


(6.9) 


Now  choose  the  curves  ( t ,  a),  t  constant,  to  be  geodesics 
with  a  an  affine  parameter,  so  that 


5 1 y ( t , a ) / 6a  =  0 


(6.10) 


giving 


1  lil  =  111  « 1 A  + 

2  3a2  <5t  5 1 


_  y  X  ,  v  £ 
RvpA$v  v  1  1 
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(6.11) 


Using  (6.4)  and  (6.10)  gives 


=  2  TT  +  RvpAC!K:vPvAlVl5lK  +  R 


6a 


=  2  61y  1VR 


>  A  i  K 


st  -  'vpA£v''1'’  +  RvuA5|<vUvAlVl5,K  +  2R  .  -  1V1C  61(mvA) 


vyA£ 


6  t 


Since  RVy 1 V  1  ^  is  symmetric,  this  is 


1  3  3t|j  =  4R  rlvl5<$lwvA  +  R  ,,vpvplvl°lA  .  (6.12) 

2F  upA?  ft  'JPPC'lA 


For  simplicity,  neglect  squares  and  derivatives  of  the 
curvature.  Then 


1 

2 


do 


4 


^vypal  1 


aAl\lP 

6  t  6  t 


(6.13) 


and  higher  derivatives  of  ip  are  all  0(R2,  VR).  Let  la(t) 
denote  ly(t,  a= 0)  and  define  the  "geodesic  connecting 
vector"  (from  Z  to  z)  to  be  a  =  al  .  a  is  independent  of 
choice  of  affine  parameter  a.  In  particular,  if  a  were 
chosen  to  be  the  geodesic  distance  along  Zz  then  la  would  be 
the  unit  tangent  vector  at  Z  to  Zz  and  o  =  (aaaa)2. 
Collecting  together  (6.7),  (6.8),  (6.11),  (6.12)  and  (6.13) 
gives 


iM  t  ,a) 


. 
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d  a 

a8v6a  a 


6 ( vevY 


Now  consider  the  expansion  of 


$(a) 


A^(  z(  t ,a) ) vA( t  ,a)  d t 


(6.15) 


Differentiating  gives  d$/do  = 


Ap  I  v1Vvy  +  Ay61y/6t)dt  = 


f(d(fly1V) 


dt  +  (Au|v  -  flv|p)lVvW>dt 


i  .  e  . 


-  dO/da  = 


u  v 

F  v  1  dt 

yv 


One  may  show  (Appendix  9)  that  continued  differentiation 
gives ,  for  n  =  1 ,  2 , . . , 


n  n 

-  d  $/d a 


'  |Flluh,  +  (n.1)F  61dlvlM“-2) 

{  yv|A(n-l)  y v | A (n -2 )  ^ 


4.  1/  \  n  d  p ,  a ,  v ,  A ,  A  (n-3) 

+  2(n-1)(n-2)Fyv|A(n-3)RA  pav  1  1  1  1_ 


+  I(n-l)  (n-2)  (n-3)F  ,,R.Uo„  6  1  P  1 0  1 V  1  A  1  ~(  n~  4  > ' 

6  y  V  A  (n-4 )  A  pa  -rr 

6  t 


d  t 


+  0 ( R 2  ,  VR) 


(6.16) 


In  terms  of  aa  =  ala,  (6.16)  gives 


$(a) 


oo 


i 

n=  0 


n  n 

a  d  $ 


n  ! 


IT 


da 


a=0 


4>  ( t ,  a  )  d  t 


^ere  we  ignore  / d ( Av 1 v )  since  it  does  not  contribute  anything  to  61 
when  the  world-line  is  varied. 


' 
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where 


<|»(t,0)  =  A  Va  -  l  F  i  aX(n) 
a  nf0  a3|^(n) 


1  ,,a  3 

V  a  + 


(n+1) ! 


+  (n+i )  ,6aaa3  + 


(n+2 )  !  6  t 


1  p  ^  3  £  6  /  i  y  i 

2ZTRe  y6a  a  a  (rrrv  +  ^ 


6a 


Y,  1 


n+3  3 ( n+4 )  6t  ^ J 


+  0( R2  ,  VR) 


(6.17) 


This  completes  the  expansion  of  (6.1).  The  action  integral 

Oi 

has  been  expressed  in  terms  of  tensors  at  Z  : 


I  = 


L(V 


a 


a3 


a 


a 


6  a 

6  t 


a 


Fa3 | X(n) 


R 


a 


3y6 


)  d  t 


(6.18) 


where 


L(t)  =  -  m(-\p(t  .a))'2  +  e  ( t ,  a )  (6.19) 

and  the  expanded  forms  of  ip  and  are  given  in  (6.14)  and 
(6.17). 


6.3  Equations  of  Motion  in  Given  External  Fields 

The  previous  section  concerned  itself  with  the  expan- 

ct 

sion  of  the  specific  action  (6.1)  in  powers  of  Q  .  The 
expansion  may  be  carried  out  quite  generally  for  any  action 
integral  representing  a  particle  with  world-line  zR(t)  in 
given  external  fields  <f>A. 

An  action  integral 
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I 


Lo(A,  9UV.  *A(z))  dt 


(6.20) 


may  be  writ  ten  as 


I 


L(Va,  ga6,  a“.  aa,  fA(Z))  dt  +  0(R2,  VR)  (6.21) 


where  =  (<f>Ai  ,  .  ,  Ra  )  and  oa  =  6aa/6t. 

£>  A  |  a  (n )  8yo 

The  equations  of  motion  are  obtained  from  variation  of 
zM(t)  with  fixed  endpoints.  Also,  since  Z  (t)  is  any  chosen 

reference  world- line,  I  must  be  invariant  under  arbitrary 

a 

variation  of  Z  (t).  This  implies  that  by  simultaneous  varia- 
tion  in  both  z  (t)  and  Z  (t)  it  is  possible  to  derive  the 
equations  of  motion  in  various  equivalent  forms. 

Consider  a  one  parameter  family  of  infinitesimal  dis¬ 
placements  in  both  zy(t)  and  Za(t),  zy(t,e)  =  zy ( t ) +eny ( t ) , 

Oi  CL  01 

Z  (t,£)  =  Z  (t)+e?  (t),  with  fixed  endpoints.  We  note  that 
since  a  and  So  /St  are  both  two-point  vector  fields1, 
differentiation  must  be  used  with  caution,  6aa/6t  and  a  06 1  3  V  ^ 
are  not  the  same:  we  have 


a  .,6  a  A 


ja  =  6aa/6t  =  6a  ( Z(t) ,z(t) )/6t  =  a  |gV  +  a  uv  .  (6.22) 


A 


The  accompanying  absolute  variations  are 


§la  =  i_ 

6  £  6  t 


a 

'  9  Z  ' 
0£ 


.  "A  =  ?.i  3Z“ 

*  -  A  a—  5 

6  £  8£ 


1aa  transforms  as  a  vector  at  Z  and  as  a  scalar  at  z.  For  a 
discussion  of  two-point  tensor  fields  (also  called  bi tensor 
fields )  see  [ 30 ] . 


' 
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(6.23) 


=  oViz6  + 

'B3e 


6e 


Uoa' 

=  6 

|%  aa 

6a 

6t 

6e 

+  aBR  “  „VY3Z5 

6  yS  iT 


Extremizing  1(e)  =  |J2Ldt  subject  to  fixed  zv ( t  . ,  e  )  = 

L  1  1 

zu(t  ,  0),  Z 06  ( t  . ,  e)  =  Za(t.,  0),  (i  =  1,  2),  yields  for 


arbitrary  variations  3zy/8e,  aZ^/ae, 


,a 


dj_ 
d  e 


,^2 


it, 


a  l  6 


ava  6t 


az 


on 


a  e 


+  iL5aa  +  1L 


„  a - 

a  cr  5  e 


3 aa  6  t 


r.  oh 
6a 

6e 


^-aBR  “  ,VY'»ZS  +  li  ®  ii 

3a“  B  y6  3l  3%  ^1“  3e 


a 


dt  =  0 


By  integration  by  parts,  the  action  principle  then  gives 


C  -  - 

d  e 


(6.24) 


2 

r  6  L 

r  a 

6a 

r  6 

'  a  L  ■ 

^  3L  o6R  y 

V5  +  It  . 

> 

3  Zal 

k 

*  a 
[6a 

-  + 

6e 

6  t 

\ 

aVa 

+  • y  36a 

o  0 

3Y A  * 

a 

> 

yr 

CX 

for  arbitrary  aZ  /a£,  az  /ae,  where 


d  t 


6  L 


a  L 


~  a  .  a 

6  a  a  a 


6  t  ^  a 


aL 


a 


(6.25) 


and  sa  /5e  is  given  in  (6.23). 

Choosing  aZa/a£  =  0  (reference  world-line  Za(t)  held 
fixed,  actual  world-line  zu(t)  varied)  therefore  gives,  from 
(6.24),  the  equations  of  motion  in  the  form 


6_L 

6  a 


a 


0 


(6.26) 


' 
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The  equations  of  motion  may  also  be  obtained  in  an  alternate 
form  by  simultaneous  variation  in  both  zy  ( t )  and  Za ( t )  , 
subject  to  <5 aa / <5 e  =  0.  From  (6.24)  this  immediately  gives 
the  transl at ional  equations  of  motion 

6  P  /  6 1  =  .  S3yV6  +  M-'i',  (6.27) 

a  2  3  y  o  a  A | a 

where 


P 

a 


3_L 

3Va 


mA  _  3  L 
M  ~  =  — 


A 


(6.28) 


Also  from  (6.24),  we  see  that  choosing  8zy/3e  =  0  (zy(t) 

CX 

held  fixed,  Z  (t)  varied)  implies  that  L  satisfies  the 
ident i ty 


(6L/6a3)a3,  -  6P  /St  +  ,  SByV5  +  M~ Y  ,  =  0 

v  I  a  a  2  3y  6a  A | a 


which  confirms  the  equivalence  of  (6.26)  and  (6.27). 

To  obtain  the  spin  equations  of  motion,  note  from 
(2.11)  the  following  identity  satisfied  by  L: 


3L  -  2—^7  g  +  ik-  +  kL_  aP  +  ~)  ^ 

,a  Z3g  Jay  _  a  „*a  3^.  A  a  B 


=  0 


(6.29) 


3  V 


3y 


3a 


3a 


which  implies 


s 

§T 


3  3  L 


a 


=  Ik  V3  +  2 

„  t,a 


3  V 


3L 


3  9 


6L 


Sy 


ay 


6  a 


a 


—  ( i  ?) 

3  ^  v  A  ;  a  B 
A  ~ 
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Antisymmetrization  and  (6.28),  (6.26)  then  give  the  spin 
equations  of  motion 

^-5 S06 3 / 6 1  =  P[aVB]  +  M~  (i  ?)  [a3]'FT.  .  (6.30) 

2  A  B 

The  derivation  of  (6.27)  and  (6.30)  closely  parallels 
that  of  (2.23)  and  (2.22).  However,  no  introduction  of  spin 
coordinates  e^  (t)  and  their  subsequent  variation  was 
needed  in  this  section;  the  spin  equations  (6.30)  are 
essentially  a  consequence  of  the  translational  equations  in 
the  form  (6.26).  Appendix  10  formulates  the  present  section 
in  terms  of  a  tetrad  field  and  shows  that  (2.22)  and  (2.23) 
follow  from  the  more  general  considerations  of  this  section. 


6.4  Gravitational  Field  Equations 

The  curved-spacetime  generalization  of  (5.48)  may  be 
obtained  by  adding  free-field  Lagrangians  for  the  electro¬ 
magnetic  and  gravitational  fields  to  the  expanded  form 
(6.19)  of  the  particle  Lagrangian  and  then  calculating  the 
gravitational  field  equations  from  variation  of  the  metric. 
This  results  in  the  one-particle  limit1  of  the  equations  of 
Chapter  4  which  described  a  continuous  matter  distribution. 


!We  may  write  I  of  (6.21)  as  I  =  fL1(x)d4x  where  Lx(x)  = 

JL ( x ) 64 ( x-Z ( x ) )dx  (x  is  proper  time  along  Za).  If  Ua(x), 
aa(x)  are  any  smooth  functions  of  x  reducing  to  Ua(x),  aa(x) 
on  Za,  and  if  Ln(x)  denotes  L(Ua(x) ,gae(x) ,...),  then  L^(x) 

=  /TjTn(x)  Ld(x)  where  number  density  /rjn(x)  =  |54(x-Z(x))dx. 
This  is  the  one  particle  limit  of  "dustlike"  matter.  (We 
consider  s4(x-Z)  to  be  a  scalar  density  at  x  and  a  scalar  at 
Z.  ) 
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The  gravi tational 
given  by  (4.26), 


field  equations  are  (4.24)  with  t  a,  $pTCT 

p 

(4.36)  where  t“B  .  and  $“By  have  the 

(mat)  (mat) 


form 


f.a8 

(mat) 


Pa(r)UB(x)  64(X-Z(t))  dT 


pa  By 

(mat) 


SaB(x)U6(T)  6  4 ( x-Z ( t ) )  dT 


(6.31  ) 


6.5  Einstein-Lorentz  Theory 

For  a  multipole  particle  in  an  E i nstei n-Maxwel 1  field, 
equations  (2.25)  to  (2.28)  inserted  into  (6.27)  and  (6.30) 
(i.e.  (2.23)  and  (2.22))  give  the  equations  of  motion 
(2.29),  (2.30).  Specific  expressions  for  pa,  SaB, 

and  q  ByS,  that  generalize  (5.24),  (5.30)  and  (5.32)  are 

a 

given  simply  by  differentiation  of  the  Lagrangian  (6.19). 
Noting  that  (-iM^  =  ds/dt,  equations  (6.14),  (6.17),  (6.19) 
give  the  following: 


a 


m 


d  t 


d  s 


,,a  ,  •  a  ,  D  u 
V  +  a  +  K  a 

Yp  o 


I 

n=  0 


3  L  /  3  V  - 

a 

rc> 

p 

ii 

«  (  VB+fo6)' 

CO 

V  e  ta 

3  Y  (n ) 

Ln  (n+1) !  6 

n  =  0 

Y  (n) 

n  yo.  8  e  6  y  (n) 

D  1  /  \ R  a  a  a  a~ 

Y3 1 Y (n)  e  6 

,  (6.32) 
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Sa3  =  2a[a8L/8a 


6] 


2m^a[a(V3]  +  a6]  +  £  6  a  £a  6  (  2  V  Y +a  Y)  ) 


V  2  e ( n  +  1 )  _[ar3]  6  y(n) 

i  (n+2 ) !  °  F  6 | y ( n ) °  a' 


n  =  0 


n=  0 


e  ar?  n  3  ]  n  e  6  y(n)  £  QQ 

T7 — 7~T\ — r  a  F  |  ,  J  ,aaaa~  c.oo; 

3  (n  +  4 )  n  !  i]  y  n  5 


m 


“6*(n)  =  23L/3F 


a3 I y (n) 


2e  _[a/.,B]  (n+1) • 3 ] \  Y (n) 

(n+l) !  0  (V  +  (n+2)  .  )a' 


(6.34) 


n ! (n+3 )  e  y 6 


R[a  .aB1aeo6(VY  +  ^Y+'aT)o 


(n+3) i y \ „Y (n) 
3 ( n+4 ) 


and 


qaBYS  =  3  L/ 3  R 


a  3  Y  5 


(6.35) 


lmdt  atafv6]v(Y  +  |oB]v[y  +  +  ioB]a'^aB1 

2  '  d  S  ^  3  J  J 


e  v  1  la  p  (  6  1  (  _  aea;f(n) 


-  f  I 


n=  0 


n  : 


!  e I Y  ( n) 


1  \i)+ 


(n+3)  3 ( n+4 ) 


(pa  and  Sae  are  parameter- independent ,  m0^-^  and  qa6y6  are 


■ 
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not,  set  t  =  x  to  fix  them  uniquely.) 

(6.14)  and  (6.17)  are  expanded  only  to  "first"  order  in 
Ra^g,  i.e.  neglecting  squares  and  derivatives  of  curvature. 
It  follows  that  (6.32),  (6.33)  and  (6.34)  are  to  first  order 
in  R  gy(S  while  (6.35)  is  to  zero  order .  An  expansion  of  I  to 
all  orders  of  R  would  enable  t0  pe  calculated  for 

any  n . 

It  was  noted  in  the  previous  section  how  the  addition 
of  free-field  Lagrangians  for  <^A  and  gag  to  a  Lagrangian  of 
the  form  (6.21)  gives,  on  variation  of  gag,  the  gravita¬ 
tional  field  equations  of  Chapter  4  with  the  material 
tensors  having  the  one  particle  forms  (6.31).  Setting  <j>  = 

(A a,  Fag),  the  whole  of  Chapter  4,  Section  4  applies,  with 


y/-g 


/-"g  J  (  x) 


eUa  64 ( x-Z ( t) )  dx 


Ta3 
(  ma  t ) 


(X) 


y 

(mat) 


(X) 


m 


aP'Y  (n) 


(X) 


■  !o(nl1 


a3y(n) 


(T) ) 


~g 


a3y6t (n) 


I  o  (  ella  ) 

=  I Q ( Sa  3  UY ) 


(x) 


agy6e(n) 


(6.36) 

(t)  )  , 


where  pa ,  Sap ,  etc.  are  the  single  particle  four -momentum, 
spin  etc.  Specific  expressions  (6.32)  to  (6.35)  may  be 
substituted  into  (4.47)  to  give  the  curved-space  generaliz¬ 
ation  of  localization  (5.48). 

For  a  collection  of  particles,  each  moving  as  a  test 
body  in  the  overall  electromagnetic  and  gravitational  fields 
without  direct  particle  interaction,  a  sum  of  actions  of  the 


. 
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form  (6.21)  will  give  the  gravitational  field  equations 
(4.47)  with  Ja  ,  T^at)  etc.  given  by  summation  of  the  right 
hand  sides  of  (6.36).  (4.47)  therefore  apply  to  a  gaseous 
assembly  of  collisionless,  charged,  spinning  particles.  In 
terms  of  an  invariant  distribution  function  y(xFn)  equations 
(6.36)  are  replaced  by  phase  space  integrals,  for  example 
<Ja(x)  =  Jy  (  x  ,n  )e(ft  )Ua  (n  )dfi  1  .  In  (3.38)  ua  is  a  timelike  left 
eigenvector  of  t^^  and  the  material  spin  flux  = 

SaSuy  is  convective,  "carried  along"  by  ua .  For  a  gas  the 
spin  flux  =  /ySapUYdn  will  not,  in  general,  be 

convect i ve . 

6.6  Propagators  for  Momentum  and  Spin 

To  compare  with  the  definitions  for  momentum  and  spin 
given  by  Dixon  [11],  in  particular  (5.1)  and  (5.2)  of  [11c], 
we  now  show  that  pa  and  SaB  may  be  related  to  pX  =  muX  via 
two  propagators  (two-point  tensor  fields  KXa,  HXa). 

Let  a-1X  (Z,z)  denote  the  inverse  of  aa,  : 


a 


3  -IX 

a  i  ,  o 


(6.37) 


a 


a 


Then 


•  a 

a 


a 


a 


(6.38) 


gives 


(6.39) 


a 


a 


where 


iCf.  [6c]  which  postulates,  on  the  basis  of  the  balance 
laws,  the  field  equations  for  a  dipolar  gas. 


. 
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„X  _  -IX  3 

N  =  -O  „CT 


a 


3  a 


HX  =  a'1X 
a  a 


(6.40) 


Substitution  of  (6.40)  into  (6.7)  gives 


iKt,a)  =  Q,  (KX  Va  +  HX  aaw,/y  ''3 
“Xy  a 


) ( Kw0  V6  +  Hy.oB)  .  (6.41) 

a  3  p 


Since  Kx  (Z,z)  and  Hx  (Z,z)  are  determined  solely  by  the 

a  a 

gravitational  field  and  the  points  Z  and  z  we  have  from 
(6.41 ) 


3  tJj  /  3  Va  =  2  q  vyKX  ,  3  ^  /  3  a“ 

r  “Xy  a 

so  that 

3  (-m(-^)^)/3Vct  =  m( -^) '^9Xy  vuKXa 

and 

9  ( "2) /3aa  =  m(-ip) '^gx^vyHXa 


2  g  vyHX 

aXu  a 

KX  p,  (6.42) 

a  X 

Hx  p,  (6.43) 

a  X 


We  now  find  expressions  for  3 4> / 3 V a  and  3<j)/3aa  in  terms 
of  Kx  and  Hx  .  From  (6.16)  for  n  =  1:  ~d$/da  =  fFx  vxlydt, 

a  a  M 

and  from  (6.15)  we  have 


V -o 


$  (  a )  =  $  (  0)  + 


d  $ 
3V 


d  v  = 


v  =  0 


• 

v  = 

» 

A  Va  dt  - 

d  t 

a 

✓  • 

J 

F_  (  z  ( t  ,v)  )  vA  ( t  ,v)  1  y  ( t  ,v)  dv 


v  =  0 


where  v  is  an  affine  parameter  along  the  geodesics  Zz. 
Inserting  (6.39)  into  the  above  yields 


I 
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$  (ct)  = 


<j>  ( t  ,ct)  dt  = 


( Va A  +  Va0  +  aan  )  dt 

a  a  a 


with  0  and  n  defined  as 

Ut  Ot 


v  =  a 


0a(t)  3  - 


FX|j(z(t,v)  )KXa(t,v)ly(t,v)  dv  , 


v  =  0 

V  =  CT 


n  (t)  e  - 

a 


F  (z(t,v))HX  (t,v)ly(t,v)  dv  . 

Ay  a 


v  =  0 


(6.44)  immediately  gives 

e  3  4)  /  3  Va  =  eA  +  e©  ,  e3cf>/3aa  =  en 

a  a  a 

From  L(t)  =  -m( -ij, )  %+e<f,  and  (6.42),  (6.43),  (6.46),  we 

p  and  Sa^  in  terms  of  and  HX  : 

a  a 

p  =  3 L / 3 Va  -  eA  =  Kx  +  e©  , 

a  a  a  A  a 

SaB  =  2a[a3L/3;  .  =  2a[a(H  B1pX  +  enB])  . 

P  J  A 


(6.44) 

(6.45) 

(6.46) 

obta i n 

(6.47) 

(6.48) 


(6.47)  and  (6.48)  are  the  definitions  of  four -momentum  and 
spin  proposed  by  Dixon  ((5.1)  and  (5.2)  of  [11c]). 
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Appendix  1 


For  a  relative  tensor  $ 

A 


U  ) 

a  a 


a 


m+1 


a 

m 


.  a 

n 


we  have 


(  1  .  ~  )  ' 
A  a 


m  10-,  3  .  3  .  Q  3  .  , 3  a  a 

V  1  r  i-l,  r.  1o3  i+l  *  m  m+1  .  n 

/  00  •••0  v.  0  0  ♦••0  0  •  •  •  o  — 

a  a  -  a  .  _  a  a.  a  .  ,  a  0  . ,  3 

i=l  1  l-l  i  l+l  m  m+1  n 


n  ,  3 ,  3  a  a  .  n  a  .  a  a 

+  l  6 b  6  1...6  m6.m  1  1(«  +  B  ) 5  p 1  ...S 

i=m+l  3  “l  %  Bm+1  Bi-1  “  Bi  61+1 


n 


n 


b/i 
-  w<5  5 

a  a  ^ 


3  a  a 

„  m  „  m  +  1  .  n  .  0 

a  0.1  0a 

m  m+1  n 


When  $A  =  (Aa,  F  Br  |  x  (n) ,  R^yfilUn)  ;  n  =  0,1,...)  we  have 


o 

MA/T  Bn  0a  3  L  ,  ~  1.0  . 

M~  ( I  ~  )  =  va  (“6  5  )AQ  + 

A  a  B  3  A  a  a ..  6, 

a  11 

00  n  ,  a  t  .  .  a  0  -  3  .  3  .  Q  3  .  ,  ,  3  „ 

v  v  1  1  n  ^  1  ri  —  l/rirPvp1i+l  o  n  r 

)  )  w™  3  **3  (-0  0  )  3  .  . 3  h  | 

^  i  2  a,  a  .  _  a  a.  a  .  ,  -  a  0  0  0.  .  0 

n  =  2  i=  1  1  i-1  l  l  +  l  n  1  2  1  3  n 


+  I  q 


a  ....  a  .  rn-1 
1  n  - 1 


n 


a 

=  4  n 


3-,  0  .  ,  0.o  3  .  0  ,  a 

l  6  1 ...  6  1  1 ( -6  166  )6  1  .  .  .6  n  n 

i  a ..  a  .  .  a  a.  a  .  .  a  0 

i=l  1  l-l  l  l+l  n-1  n 


3 1  3  a 

+  6  1.  .6  n  1(S  n6p  ) 

an  a  ..  a  0 

1  n-1  n 


n 


0.  0.0.  0  .  .0 

123  4  n-1 
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oo 

1  y 

n 

y 

a 

2  ^ 
n  =  2 

L 

1=1 

n-1 

„  “l 

.  .  a 

Cl  -I  •  •  Cl  #  M  Cl  #  -  •  •  Cl 

m  1  1-1  1+1  n  r- 

m  F 


n  =  4  i=l 


'a 


.a  .a 
n  - 1 R  n 


n 


a,  ...  a  .  - aa  .  ,  -  .  .  a 

1  1  l-l  l+l  n 


n  -1 


°^..|  •a-^_xaa-j_+x*  •  •  a 


00  o.  i  .  .  .a  -  0 

1  n-1^3 

,  'a  a_  a. a_  a ....  a  - 

n  =  4  1  2  3 1  4  n-1 


+  l  q, 


The  symmetries  of  q  and  m  together  with  pa  =  Pa 


M~(I  ~)  =  (p  -  P  )u3 

A  a  B  a  a 


l  mB^(n)F 

n=  0 


,.M  -  y  I  ( n  +  1)  mYy  ^  (n)  ®  p 
ay|X(n)  2n=0 


I  q 


yy  v  X  (n)  ^3 


n=  0 


a 


yy v | X (n) 


3  l  q  ^U")r 


n=  0 


00  /  \ 
l  (n+l)  q  YWA(n)SRS 


n=  0 


yy  v  |  ( X ( n) a) 


-eAa  give 


yy I (  X (n) a) 
5 

ay y ] X (n) 


which  inserted  into  (2.22)  gives  (2.30). 
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Appendix  2 

Consider  a  class  of  Lagrangians 


L  ( u  ,  3  a 

a 


m 


,  e  )  >  l(u“. 


3  am , 

a 


a 

,  -u  u  , 

a 


a .  m  \ 
u  3  a  ) 

a 


where  L  is  an  arbitrary  function  of  its  last  two  arguments, 
subject  to  L ( ua , 3aam , 0 A , 1 , 0 )  a  fixed  function  of  ua,  0A  and 
3oam.  (e^  is  shorthand  for  (a"\e<a)  ,e^  .  *A'*a|«  ■  R%Ys  >•  >  Let 

AX  =  6X+u  uX(-u  uT)  1.  From  (any)  L  we  can  form 

a  a  a  y  J 


L'(ua,  3  am ,  e  )  =  L(ua/(-u  uX)^2,  Ax3  am,  e  ) 

a  A  A  a  A  A 


=  L(ua/(-uxuX)^,  AX  3  x  am  ,  0  A ,  1,  0) 


so  that  each  L  gives  the  same  L'  . 

We  will  show  that  variation  of  the  metric  for  any  L 
gives  a  total  energy  tensor  T  a  that  is  expressible  in  terms 

of  L'  .  Since  each  L  determines  the  same  L'  it  then  follows 

a 

that  Tp  is  independent  of  the  particular  choice  of  L. 

To  calculate  Tpa  from  L  we  follow  the  same  steps  as  in 
sections  3  and  4,  modifying  (3.11)  to  include  ua. 


I  = 


L  <  ^  A  ’  ^  A  I  a  ^  d4><  +  (16,T) 


-1 


/-g  R  d4  x 


,  _  /  a  m  _  (a)  ,  na  \ 

^ A  (u  ,  a  ,  ea  ,  4>a,  R  ^6) 


The  variation  in  ua  resulting  from  variation  of  the  metric 
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may  be  computed  from  (3.1)  and  (3.3).  Since  nya3y  = 

(-g)  2eya3y  where  eya3y  is  the  Levi-Civita  permutation 
symbol,  (3.3)  implies  that  6Ny  is  parallel  to  Ny ,  implying 
<$uy  parallel  to  uy  ,  <$uy  =  euy  say.  6(gaeuau3)  =  0  then 
implies  <$gagu  u  +2egagu  u  =  0 ,  so  that  2e  =  6ga^u  u  and  suM 
=  y5ga^u  u  u  .  This  gives  an  additional  term  ( 3L/ 3uy ) uyupua 
in  the  r.h.s.  of  (3.20).  With  ^ A  containing  ua  the  funda¬ 
mental  identity  gives  an  additional  (3L/3up)ua  on  the  r.h.s. 
of  (3.22).  The  above  two  remarks  imply  that  the  field 
equations  (3.24)  contain  an  additional  ( 3L/ 3uX ) ApUa  on  their 
r.h.s.  so  that  (3.26)  and  (3.27)  are  replaced  by 


t  0  +  ( 3l_/3uX)  AXua  =  t  0  .  +  t  ° 

P  P  p(mat)  p(40 


X  ..o  ,  „  _mN  „  i  _  ( a) 


( A2 . 1 ) 


p  (mat) 


L.fi°  +  3L  Au°  -  (3  am)  3  L 

1  n  - o  o  - 


-e 


3  U 


p  '  - —  a  |  P  - 

3  ( 3  am)  3  e 

o  ala 


3L 

TO 


and  tp^  given  by  (3.28) 


a 


To  express  t  "  in  terms  of  L'  we  note1 

p  (  mat) 


3L'  3  L  3 

a 


3  U 


3U3  3  U  a 


8 


7  TO 

~  U  A  U  )  2J 


3L 


3 ( 3 R am)  3ua  3  X 

p 


—  (<5  3  +  u  U3)  + 

3  u  3  '  a 


3  L 


(  X  ^ 

u  u 

e 


3 ( 3  0  am)  3  u  a 

p 


-u  u 

Y 


Y 


3 .  a 

A 


m 


’We  are  not  concerned  about  the  arbitrariness  of  tpa  or  Tpa 
as  functions  of  ua  and  3aam  due  to  the  constraints.  (-uau^) 
and  (ua3aam)  are  therefore  set  equal  to  1  and  0  in  the 
following  formulae  after  differentiation. 


■ 

' 


■ 
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-  *V3 

.  p  a 

a  u 


+ 


— — -  (  6XU  + 

a(aQam)  a  3 
p 


uX g  +  2  u  uXuQ ) 3  a 

w  3a  a  3  X 


m 


i  .  e . 


a_L_'  =  aL  a3  + 
aua  a a  a(a0a“l) 


u  ( a  am) 
mN  3  a 


(A2.2) 


3 


Simi 1 ar 1 y 


aL' 


3L 


a ( a  am)  a  (a  a11)  a  (a  am) 

o  y  a 


(A.Van) 


sJ. _  A°  .  (A2.3) 


m\  y  x 


a ( a  am)  Y 
y 


According  to  (A2.2),  ( A2 . 3 ) ,  L'  satisfies 


3L#  a  n 
-  u  =0 

a  ua 


3L' 


3(3  am)  a 
a 


u  =  0 


( A2.4) 


(A2.2)  and  ( A2 . 3  )  inply 


ULaV  -  O  am) — it 

Xp  p  „  ,  „  _  m 


a  u 


a  (a  a  ) 

a 


AX  +  u.O  am) 

a  ux  p  a ( a  _  am)  3  p 

3 


u°  , 


- ( 3  am)_  3L_  a°  =  xtl  u°  -  (9  am)  3‘-' 

~  ~  p  ,w  *  „m 


p  ’  a  ( a 0 am)  3  aup 
3 


a  (a  a  ) 

a 


(cf.  (3.32)).  These,  together  with  (3.27)  and  (3.32),  imply 
that  tp^mat)  of  (3.27)  and  (A2.1)  are  the  same. 
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Appendix  3 


rv 

The  contribution  to  P  from  L ,  .  N  =  enA  u  /(-u,u  )2  is 

p  (A)  a  X 

P  ...  =  aL,../aup  -  L/aNu  =  enA  +  enA  uau  -  L,.,u 

P  ( A)  (A)  (A)  p  p  a  p  (A)  p 

enA  . 


The  contribution  to  P  0  from  L 


(A)  18 


P  ? 4i  =  Lm4°  '  L,..n  1(3  am )  ( a n/ a  ( a  am  )  ) . 

p  (,  A)  ( A)  p  (A)  p  a 

Formulae  ymn  =  gyv(a  am ) ( a  an ) ,  ymny  =  5m,  a  = 

y  v  np  p 

Y„„ ( 3„ am )( 9 v an ) ,  y  =  det(ymn)  and  n  =  Ny^  (cf.  section  2) 


mn  y 

imply 


( 3  alu)(an/8(a  am)  )  =  N(a  au‘)(yy  2)(yy  )a(yr  ) /a  ( 3  a“l) 
p  a  p  l  Pn  a 

H  m,  a\)„  n  a 

=  Ny  Ymn(3pa  )g  3va  =  nAp. 


m  w  1 


pn 


m . 


Therefore  P  ,  .  *  -  0 . 

p  (  A) 
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Appendix  4 


6(e)L=  >-AS(n)«(e)(*A|o(n)) 

n=  (J  ~ 


-  L  5(e)V  +  l  L' 


Aa ( n) x 


n=  0 


^(e)^A|a(n)x^ 


( A4 . 1 ) 


Equation  (2.8)  implies 


.  Aa (n) x  ,  s 

L  ~  8/  \  /  \  ) 

(e)  A  a  (  n)  x 


Aa ( n ) t 


n 1  .  / T  B§ (n) .  a  *rP  ■) 

(e)  vyA|a(n)  I  t  A  a ( n )  p  ^  B  |  (3  (  n  )  a  x  - 


({5r0^ (V 


=  _  |_A~(n)T,  6  (xl>  |  )  +  (div) 

x  (e)  A  a (n) 


Aa(n) t  ,  B N  a 


,  Acra  ( n-1)  x  ^  p 

A  'p  rB[a(n)  ^A | (pa (n-1) ) '  ax 


according  to  repeated  application  of  (2.5).  If  we  repeat  the 
above  procedure  n  times  we  obtain 


.  Aa (n) x  ,  .  v 

L  ~  o,  ,  ) 

(e)  A  a  ( n  )  x 


n  +  1LAa(n)T  ,  *  +  (dlv) 

x  a  (  n)  (  e)  A 


V  /  t  N  m  f,  Aa  (n-m)  3  (m)  x 

+  L  (-1)  iL  ~ 

m=  0 


/  x  a  , 

3(m;  a.  P  B|a(n-m) 


,  N  .  Aaa ( n-m-1) B (m) x 

-  *  '  iBU^Al  (pa(n-m-l)))6r^T 


Inserting  this  into  (A4.1)  and  noting  that 
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oo 


l 

n=  0 


n 


I 

m=  0 


00  00 

I  I 

m=0  n=m 


OO  00 


l  l  ( r=n-m) 
m=0  r=0 


gives 


6 


(e) 


L  = 


l  (-l)nLA~ ^ 

n=  0 


a  (  n)  ^  (  e)  ^  A 


+ 


(  d  iv) 


OO  00 


+  I  l  (-1)  (L 

m=0  r=0 


m  r ,  Aa ( r ) 3 ( m)  x 


(i  B'i  a  i 

3(m)  A  p  PB  |  a ( r) 


-  r 


Aaa ( r-1) $ (m) t 


* 


3  (m)  A  |  ( p a ( r -1 ) ) J  ax 


1  5T1 

*  i 


=  L*S(e)'(,A  +  (div) 


+ 


00 

I  (Lt?(r)T(IAB) 

r  =  0 


a(r) 


■ Affa(r-l) t , 
rL*  ' 


( pa ( r-1) ) 


P 

a  x 


=  LA6,  *ip  +  UTa  6TP  +  (div) 
*  ( e)  A  pax 


according  to  (4.3),  (4.11). 
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Appendix  5 


Expand  f(z)dzp/dx  in  powers  of  aa: 


f(Z+a)dz3/dx  =  (UB  +  a3)  I  ~X(n) 

Lj  -n  ' 


n=  0 


ns  a~  '  3X(n)f(Z)  <A5-1) 


W  .  . 

x  x  V  1  ,„B_Y  x  w_if(Z)  . 


=  UBf(Z)  +  ap  f ( Z )  +  l 


n=  0 


WT<U  a'  + 


To  form  the  antisymmetric  part  of  the  summation  in  the  above 
equation,  consider 


d_ 

dx 


l  ,-rrrvT  aBa-(n)3,  ,_,f(Z) 


n=  0 


( n+1) ! 


X(n) 


v  1  r*3  X(n)_  r  ,  B*y  X(n-1)„ 

I  77ZTTi(a  3-  '  'f  +  no  a  a~  8 


n=0 


(n+1)  ! 


X  (n) 


■p  X(n)MY0  f) 

YX(n-l)f  +  °  °'  U  \X(„) f) 


i6f(Z)  +  I 
n=0 


1  •  3  Y  ,  (n+1)  B*y  ,  1  _ 3 1 |Y 

T^2T!a  a  +  T^W'°  0  +  ■aSiY.0  U 


-(n)V(n)f(Z) 


(A5.2) 


Using  (5.24)  and  subtracting  (A5.2)  from  (A5.1)  gives 


f (z) dz3/dx  -  d 


dx 


l  aBa~ (n) 3,  ,_N  f(Z) 


n=  0 


(n+1) ! 


X  (n) 


J 


(A5.3) 


=  f  C  z )  U 6  -  e’1  l  mBY^(n)(x)3YX(n)f(Z) 

n=  0 
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This  identity  is  given  in  [5,  page  242].  When  used  in 
[5]  to  obtain  the  translational  equations  of  motion  (equa¬ 
tion  168,  page  201)  the  total  proper-time  derivative  is  not 
combined  with  mu  ,  as  we  do  in  (5.30),  but  is  placed  on  the 
r.h.s.  of  the  equations  as  part  of  the  total  four-force. 

A  short  derivation  of  the  polarization  equations 
(5.17),  (5.18),  is  available  from  (A5.3).  Setting  f(z)  = 
es4(x-z)  then  integrating  w.r.to  t  immediately  gives 


.  a  ,  \ 

J  ( X)  = 


edz/dx  64 ( x-z(t) )  dx 


n  +  l 


eU“s4(x-Z(T))dT  -  l  (-1)  3  gx(n) 


n=  0 


X  ~ 


m 


aB-(n) (t) <$4  (x-Z(t) ) dx 


=  J“(x)  +  3BHa6(x) 


according  to  (5.25)  and  (5.26). 

If  we  set  f(z)  equal  to  64(x-z)  in  (A5.3),  then  mult¬ 
iply  the  resultant  identity  by  mua  and  integrate  w.r.to  t, 
we  obtain 


ta3(x)  = 


muadz^/dx  64(x-z(t))  dx 


muaU364(x-Z(x))dx  -  e"1  l  (-l)n  9  y^(n)J 


-1  r  ,  ..n+l, 
n=0  X 


muam^Y~('n')  (x)64(x-Z(x))dx 


l  3 

n  =  0 


(n  +  l)  !  uxHn)  J 


'  d  ( mu  a_)  g^a~(n)(54(x-Z(x))dx  . 

dx 


. 
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From  (5.20),  (5.24),  (5.22)  and  (5.23)  this  is 

_  -j-a  3  +  3  +  /\a  $ 

Y 
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Appendix  6 


Expanding  aaf(z)dzy/dT  gives 


Y 


aaf(z)dzY/dT  =  aaf(Z+a)(UY  +  aY) 


oo 

I.  TTT  f<z>,A(n)°a°-(n)(uY  + 


n=  0 


Let 


XaY  =  aaf(z)dzY/dx  -  l  e  Xf(Z) 

n=  0  +  (n) 


(A6. 1 ) 


maY-(n)(T)  .  ( AS . 2 ) 


From  (A6.1)  and  (5.24)  we  obtain  X 


ay 


l  f(Z)  Wnl0 

n=  0  ’^(n) 


X  (n) 


( 


1 


n  !  (  n  + 1 ) 


,)a“U 


a  1 1 Y 


r  1  (n+1)  ^  a*y  ,  1  y,,a  (n  +  1)  ya 

+  b-(S!)!)J  a  +  (n+TT'0  U  +  ( ^T+2T !  ° 


=  l  f(Z), 


n=0 


X(n) 


a 


X(n) 


n  a, ,y  .  1  y,,a  ,  (n+1) 2  a*y  ,  (n+1)  ya 

(^iy!a  u  +  (^!°  U  +  (n+2)~ ! °  °  +  (n+2)~ ! 0  ° 


(A6.3) 


When  f(z)  is  set  equal  to  (z)  we  find  from  (5.30)  that 

Y 

the  second  term  in  the  above  summation  is  just 
e_1(mue-pe)Ua.  This  term  will  combine  with  mu^aU6^  of  (5.31) 
to  form  p[aU^ .  We  wish  to  express  the  remainder  of  (A6.3) 
as  a  combination  of  derivatives  of  f(Z)  coupled  to  multipole 
moments  ma$Y(n),  together  with  a  derivative  w.r.to  t.  It  is 
appropriate  to  consider  the  derivative  w.r.to  x  of  a  term 

having  the  form  £  anf  ^  ( Z ) aaaYa~ ^  where  an  are 

n=  0 


' 
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numerical  coefficients.  From  (A6.3)  the  most  simple  choices 
for  an  are  ( n+ 1 ) 2/ ( n+2 ) !  or  ( n+ 1 ) / ( n+2 ) ! .  Making  the  latter 
choice  we  rind  from  (A6.3)  and  (5.24)  that 


X“Y  -  I 


1  A(n)  Y,,a  d 

/  / — rArriT(A)  w  s a ~  a  U  -  -r— 

(n+1) !  ’A(n)  dx 


oo 

V  (n+1)  a  y  A(n) 

&  (^«T!f(Z).A(n)a  0 

n=U 


uu 

£  f (Z)  (n) ~a 


n  =  0 


A(n) 


a  ~  a 


n 


1 1 Y  ,  f  (n  +  1)  2  (n  +  1)  'i  *  Y 
(n+1) ! U  l(n+2) !  (n+2) ! J 


n  =  0 


(n+1)  a  y 
(n+2) ! °  ° 


(f(Z) 


|5  X  (n) 


A (n-1) • ^ 


.  .  U u  a  ~  v  “ y  +  n  f  (  Z)  .  s  ct~  '  'a 
’  A ( n) 5  5  A ( n ) 


n' 


I  f(Z)  0“a-(n)a'5 

n=  0  +  (n)« 


(  n  +  1)  .. y  ,  (n  +  1)  (n  +  2  )  •  y 

(n+2)!  ( n+3 ) ! 


/  \ 

V  r/7N  a  A (n)  Y 

L  f(Z).X(n)6a  0 


n=  0 


(n  +  1)  .,6  (n  +  1)  (n+2)  *6 

(n+2) !  (n+3) ! 


I  e‘1(n+l)f(Z)! 

n=  0 


m 


6  Ya A ( n) 


( A6.4) 


(A6.2)  and  (A6.4)  give  the  required  expansion: 


aaf (z) dz  r /dx 


Y 


I  e_1f(Z) 

n=  0 


»  A  (n) 


m 


a  y  A  (  n  ) 


+  l  e  1(n  +  l)f(Z) 
n=  0 


»  A ( n) 6 


m 


6  Ya  A ( n) 


+ 


n=0 


(n+1) !f(Z) >X(n) 


a*(nVu"  +  d_ 


^n=0 


(n+1)  „a  y_^ (n) 

Xn+2)" !  f  Z)  >X(n)°  a 


( cf .  [ 5 ,  page  242 ] ) . 

Choosing  f(z)  =  fe  (z)  gives  an  expression  for 


Y 
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aaf ^ (z)dzy/dT  which  on  substitution  into  (5.31)  gives  the 
spin  equation  of  motion  (5.33). 


Appendix  7 


Let  AaS 


denote 


(Air)  1(fa  H3X  -  ff  fUVga3) 
X  4  yv 


The  field  equations  fag  =  4ir(Ja+Mae  ),  fr  _  , 

>3  >3  [ a  3  >  y  J 


A 


6  =  -  f  J6  -  yMXyf 

a»3  a  3  2  A  y  ’  a 


Def i ne 


00 

Ma3y(m)^x^  E  ^  (-1) nma^T (m) T ( n) 

n=  0 


»  X  (  n) 


then 


m 


a3y(m)  _  ^.a3y(m)  +  ^a3Ay(m) 


From  (A7.4)  we  have 


MXy  f  =  mXy  f  -  ( MXy T  f  )  +  MXyT  f 

Ay>a  Ay»a  Ay>a  >x  Ay  >ax 


Simi lar ly 


MXyxf. 

X  y  »  a  t 


mXy x  r 
m  f , 

X  y  j  a  x 


(MXuTYf,  )  +  HXpTYf 

X  y  »  a  x  »  y 


(A7. 1 ) 

0  imply 

(A7.2) 

(A7.3) 

(A7.4) 

.  (A7.5) 

X  y  j  a  t  y 


Repeated  application  of  (A7.4)  gives 


' 


00 


CO 


I  m 

n=  0 


Xy t (n) 


Xy  >  ax ( n) 


-  (  I  M 

n=  0 


Xy t y (n) 


^Xy ,ay(n)^ 


From  (5.35),  (A7.2)  and  (A 7. 6)  we  obtain 


3  Ta  ^ 

6  (mat) 


a 


>3 


+ 


00 

(  l  ^XyTy(n) 
n=  0 


^  Xy  >  ay ( n) ^ 


According  to  (5.37)  and  (A7.1)  this  is  just  (5.36). 
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Appendix  8 


In  the  equation  following  (5.38)  we  wish  to  express  the  two 
summation  terms  on  the  r.h.s.  as  (cf.  (5.37))  plus  a 

divergence.  From  (A7.4)  we  have 


I  fa  w  ,m 

n=  0  Y-X(n) 


B yX (n) 


I  f 

n=  0 


.a 


o°  c° 

MBYA(n)  +  g  (  l  M^(n))  -  I  f 


Y » X (n)  '  T»X(n) 


■a  j^PyT  X  (n) 


n=0 


y  5  X (n) x 


=  f“  MBY  +  3  (  I  f“  . 

y  y  ^  L n  K>X(n) 


M 


3Ky X (n) 


Y  n=0 


) 


(A8. 1 ) 


and 


^  (n  +  1) f  y  ,  5X(n) 
n=  0 


m 


y  6  B  X  ( n) 


,y  6  3  X  (  n) 


.a 


I  <^fUY.5A(n)M,UPi'1W  +  3x(.L(n+1)f,Y.«(n) 


n  =  0 


M 


y63xX(n)j 


n  =  0 


I  (  n  +  1 )  i 
n=  0 


y  >  6  X ( n) x 


M 


y  6  B  x  X  (  n) 


The  last  term  is  -  l  nfa 

L  n  y  » o  X ( n) 
n  =  U 

the  first  term  on  the  r.h.s  to  give 


which  combines  with 


l  C-i>f\,ax(n)*YSe-(n)  (A8'2) 

n  =  0 


V  ,a  MY«BX(n) 

in  Y.«i<«) 


3  (  l  ( n+1) f 

Y  n=  0 


a 


6X  (n) 


H 


k  6  B  y  X (n) 


) 


Substitution  of  (A8.1)  and  (A8.2)  into  the  equation 


■ 


■ 


following  (5.38)  and  use  of  (5.37),  (5.39)  gives  (5.40). 
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Appendix  9 


From  section  6 . 2  we  have 


and 


$  (a) 


Ax(z(t,a)  )vX(t,a)  dt 


-  d  $  /  d  a 


F  vulV  dt 

y  v 


Continued  differentiation  and  use  of  (6.4),  (6.10)  gives 
(neglecting  derivatives  of  curvature  and  squares  of 
curvature) : 


- d  2  $ / da 2 


(F 


y  v  A 


y  v  ,  A 

v  1  1 


+ 


F  1 

y  v 


v 


61 
6  t 


y 


)  d  t 


-d3$/da3 


(F 


yv| Ax A2 


vylvlXllX2 


+  2  F  ,.lVlX6ly 

yv|x  ^ 


+  F  R,U  vplalV 

yv  A  po 


-d4$/da4 


f  F  |  VylVl  A  1  A2  1  A  3  +  3F  i  v  -j  X  x  i  A  2  6  i  y 

^  y  v  |  A  ^  A  2  A  3  y  v  |  A  i  A  2 


+  3  F 


R  y  vp]a1v1A1A1  +  F  R  y  i ° i v i A 6 i P] 

A  n  n  1 1  \)  A  n  n  -  J 


y  v  I  A  i  A  pa 


y  i  0  l  v  i  ^  i  P' 


y  v  A  pa 


6  t 


The  above  expressions  are  all  of  the  form  (6.16).  We  may 
verify  (6.16)  by  induction:  differentiating  (6.16)  gives 


)  dt  , 


d  t  . 
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-dn+1t/dan+1 


(FUvU(n)VyiVl'<n)  +  <1+(n-1>)FWv|X(n-l)lVl-(n'1)fiV 
+  ((n-l)+i(n-l)(n-2))F]lv|i(n_2)  R/ pa  VP  1  0  1  V 1X  1  *  (n~2) 

+  (i(n-l)  (n-2)+|(n-l)  (n-2)  (n-3)  )  Fyv  |  ^  (n.3)  R^l  ^  ^  ^  '  (n'3)||P]  «  • 


The  coefficients  in  the  above  expression  are  n,  n(n-1)  and 
n(n-1)(n-2)  so  that  the  above  is  just  (6.16)  with  n 
replaced  by  n+1 
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Appendix  10 

Consider  a  composite  particle  consisting  of  n  point  par¬ 
ticles  with  world-lines  z^it)  (i  =  1,2,  ...n).  With  a 
common  central  world- line  Za(t),  the  Lagrangian  of  each 

particle  i  may  be  expanded  in  powers  of  connecting  vectors 
a(i)’  "^e  total  action  for  the  composite  particle  in  an 
external  field  $  is  then 


I 


L  d  t  with 


n 


L  =  I  L 


1=1 


(i) 


(  A  1  0  .  1  ) 


and 


L 


(i) 


I  is  a  functional  of 


( A 1 0 . 2  ) 


XA(t)  =  (z^  (t)  ,  1  =  1 _ n)  ( A 1 0 . 3 ) 

The  4n  translational  equations  of  motion  are  obtained  from 
di/de  =  0  for  infinitesimal  displacements  of  each  z^^(t) 
with  endpoints  held  fixed.  In  other  words,  from  demanding 
di/de  =  0  for  arbitrary  variations  3X  (t,e)/3e  subject  to 

r\. 

X  (t  , e)  =  X  (t  ,0)  (i  =  1,2).  This,  however,  gives  too  much 

A  i  A  i 

information  when  one  is  interested  in  the  overall  motion  and 
spin  of  the  composite  particle  and  not  its  detailed  internal 
dynamics.  Each  particle  i  is  governed  by  equations  (6.27), 
(6.30)  and  the  equations  for  the  composite  particle  are 


obviously  given  by  summing  over  i  in  these.  The  equations 
for  the  total  four -momentum  and  total  spin 
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P 

a 


n 


3L 


(i) 


i=l  3  V 


a 


11  r 

.a  v  Let 

'  b  =  ii  a<« 


3L 


3a 


(i) 

3] 

(i) 


n 


3L 


M' 


I 

i=l 


(i) 


3Y 


(  A 1 0 . 4  ) 


are  then  (6.27)  and  (6.30). 

(a) 

To  compare  with  Chapter  2,  Section  3,  let  ea  (t) 

a 

denote  any  orthonormal  tetrad  defined  along  Z  (t).  Two  types 
of  (constrained)  variation  in  X A( t )  will  be  discussed  below 
that  are  essentially  the  variations  considered  in  Ch.  2, 

Sec.  3.  The  action  principle  appearing  there  is  then  a 
consequence  of  the  more  general  principle  di/de  =  0  for 
arbitrary  3X  /3e  of  this  chapter.  Instead  of  independent 
variation  of  each  z^.N(t),  the  action  principle  of  Ch.  2, 

( i) 

Sec.  3,  selects  just  the  variations  needed  to  deduce  the 
composite  equations. 

Let  each  have  scalar  components  o ,  w.r.to  e^a^: 

(i)  ^  (i)a  a 


a 


(i) 


=  a 


(  i)  a 


e(a)a. 


a 


(i) 


°  r  •  \  e 
(  i)  a 


(  a)  a 


0  r  ■  \  e 
(  i)  a 


(a)a.  (  A  1 0 . 5 ) 


It  follows  from  (A10.1),  (A10.2),  (A10.5)  that  L  may  be 
written  as1 


L  =  L  ( V 


a 


(a) 


a 


e(a)  . 

a 


V 


a 


(i)a’ 


a,.,  ) 

(  i)  a 


( A 1 0 . 6  ) 


1 L  denotes  the  Lagrangian  both  as  a  function  of 
( ga(3  >a?i)  >  cr^i)  )  and  as  a  function  of  ( eia)  ,  a  (  i)  a  » o  (  i)  a  ) 

since  the  particular  meaning  will  always  be  clear  from  the 
arguments  differentiating  L. 
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From  (A10.5)  it  follows  that 


3  L 


3  e 


( a)  3 


y  3L(P  8g(i) 

i=l  So).,  34(a)B 

(  l) 


n  3  L  ,  . . 

T  — ^  sic... 

i  =  l  3 aj.,  6  (l)a 

(  l) 


giving 


(  a )  a  3  L 


3  e 


( a)  3 


n 

I 

i=  1 


a 


3  L  ,  .  x 

a  (  i) 

( i) 

30(i) 


Definitions  (2.19)  and  (A10.4)  for  spin  are  therefore  the 
same.  Furthermore,  we  have  and  M~  given  by  (A10.4)  when  L 
is  written  either  in  terms  of  the  connecting  vectors  or  in 
terms  of  the  tetrad. 

Consider  the  equations  resulting  from  di/de  =  0  for 
arbitrary  3Z  /3e  with  ea  and  a(i)a  held  fixed  by  parallel 
propagation.  This  is  exactly  the  action  principle  generating 
(2.23) 1 .  di/de  =  0  is  given  by  summing  (5.24)  over  i,  and 
6e^aVse  =  0,  <Sa(i)a/6£  =  0-  imply  =  02.  The 

resulting  equations  are  therefore  (6.27)  (i.e.  (2.23))  for 
the  composite  particle. 

(  a  ) 

Equation  (2.22)  was  obtained  by  variation  of  ea 

subject  to  Sg  0/5e  =  0,  with  3Za/3e  =  0,  3a,..  /3e  =  0.  This 

J  3a3  ( a 


^he  scalar  components  a(±)a  P^aY  no  part  in  the  derivation 
of  the  equations  of  motion  since  they  are  held  fixed  here 
and  in  the  next  type  of  variation  considered.  As  scalars 
their  infinitesimal  generators  are  zero  so  they  do  not 
contribute  to  the  invariance  identities.  No  mention  of  the 
dependence  of  L  on  scalars  such  as  ^(a)a  and  electric  charge 
e  need  therefore  be  made. 

2The  resulting  variation  in  z^  is  giyen  by 
0  =  SoJ^/Se  =  a?i)  Jo3Z3/3e  +  of  ±)  I  x  3z)  i)  /  3  e  . 

For  zero  curvature  where  =  z?i)-Z  ;  ?fi>  1 6  = 

and  h  =  Sy)  this  is  a  rigid  translation 

3Z^.)fi;If/3e  =  3ZMtfc)/3E. 
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induces  a  variation  in  XA(t)  given  by  5a(i)/6e  = 


a 


(  i)  a 

di/de  = 


6e^a^a/6e  =  CT(  ±)  |  x 3z  ( i)  /  8  e  1  »  and  from  (6.24)  gives 


I|(3L(t)/3o“  )6a“  /fiedt  =  0.  Let  6e(a)/5e  = 

L±)  (i)  (i)  a 

where  <5g  g/6e  =  0  implies  =  0.  Then  6a(i)/5e:  = 


3^  ( a) 


p  "e 
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Inserting  this  into  di/de  =  0  gives 

[a  3  ] 


a 3c.  (  a)  _  3 


the  spin  equation  of  motion  as  ^5L(  t  )/<5a  =  0. 

i 

Equation  (6.29)  and  the  equation  following  it  then  give 
(6.30)  (i.e.  (2.22))  for  the  total  spin. 

To  summarize,  extremizing  I  for  arbitrary  variations 
sz^^/ae  leads  to  equations  of  motion  (6.27),  (6.30)  for 


each  i,  while  the  constrained  variations  induced  in  z 


A 


(i) 


from  those  of  Chap.  2,  Sec.  3,  give  the  same  equations 
summed  over  i  for  the  total  four -momentum  and  spin. 


^or  zero  curvature  U  =  (z(a  -Z  Kx  =  ® A  so  that 

^z  G)  (t,€)/n  =  ( t,t)/lc  =  t )  <3-  p  which  is  a 

"rigid"  rotation.  V 


